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ABSTRACT 


The  combinatorial  geometry  allows  efficient  Monte  Carlo 
radiation  transport  calculations  of  detailed  three-dimen¬ 
sional  geometries.  Because  regions  and  media  are  formed 
by  combination  of  basic  bodies  such  as  boxes,  spheres, 
cylinders  and  others,  the  input  required  of  the  user  is  both 
relatively  simple  and  easily  modified.  The  M0RSE  code  is 
a  multigroup  neutron  and  gamma  ray  transport  Monte 
Carlo  code  that  may  solve  either  neutron,  gamma  ray,  or 
coupled  neutron-gamma  ray  problems  in  either  the  forward 
or  adjoint  mode.  M0RSE  has  a  wide  variety  of  available 
input  options,  including  splitting,  Russian  roulette,  expo¬ 
nential  transform,  energy  biasing,  importance  regions, 
albedo  surfaces,  and  the  scoring  options  available  in  the 
SAMB0  analysis  package. 

This  document  details  the  incorporation  of  a  new  version 
of  combinatorial  geometry  into  the  M0RSE  code  and  is 
meant  as  a  user's  manual.  All  new  subroutines,  including 
the  combinatorial  geometry,  are  documented  in  detail,  in¬ 
cluding  flow  charts.  Input  instructions  for  all  M0RSE 
options  and  layouts  of  all  M0RSE  commons  are  given  for 
user  reference.  The  PICTURE  code  has  been  modified  to 
work  with  the  combinatorial  geometry,  allowing  two-dimen¬ 
sional  printer  pictures  of  combinatorial  geometry  models. 
Two  M0RSE  sample  problems  and  a  PICTURE  sample 
problem  are  discussed.  Input  and  output  from  these  prob¬ 
lems  and  the  complete  code  are  available  from  the 
Radiation  Shielding  Information  Center. 
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I.  INTRODUCTION  « 

(i) 

The  M0RSE  Monte  Carlo  code  was  first  distributed'*7  in 
early  1970  with  a  revised  version^’  ^  released  in  the  fall  of  1970. 

'  (4  81 

There  have  been  several  additional  features'  '  '  added  in  the  past  year. 

\ 

This  report  describes  the  addition  of  a  version  of  the  MAGI  developed 
combinatorial  geometry'  ’  (CG)  to  the  M0RSE  code  and  also  serves  to 
document  a  version  which  is  operational  on  the  UNIVAC  1108.  The 
definition  of  variables  in  labelled  commons,  layout  of  blank  common 
and  input  data  have  been  collected  for  easy  reference.  A  description 
of  the  PICTURE  program  for  the  combinatorial  geometry  package  and 
descriptions  of  routines  required  to  solve  two  sample  problems  are 
included.  , 

l 

The  replacement  of  the  05R  geometry  module  with  the  CG 
module  required  significant  modification  of  the  CG  package.  These 
modifications  were  required  so  that  the  CG  piodule  could  be  added  to 
M0RSE  rather  than  adding  MORSE  to  the  CG  package.  Only  those 
M0RSE  routines  that  have  an  interface  with  the  geometry  package  were 
changed  and  all  options  in  M0RSE  are  available  with  the  CG  module. 

To  retain  all  M0RSE  options  it  wa^  necessary  to  write  several  new 
routines  as  part  of  the  CG  module. 

V  Features  of  M0RSE  which  are  dependent  on  the  geometry 

i 

module  include  the  ability  to  (1)  deter  mine,  both  the  media  and  an  im¬ 
portance  region  given  a  set  of  coordinates,  (2)  to  track  a  particle  through 
the  system  including  internal  voids  and  surrounded  by  an  external  void, 
i  (3)  determine  the  number  of  mean  free  paths  between  any  two  sets  of  co¬ 
ordinates  (one  may  be  in  either  an  internal  or  an  external  void),  (4) 
determine  the  vector  normal  to  a  surface  at  any  point,  and  (5)  permit 
particle  "'eflection  at  a  boundary  without  "getting  lost.  ”  Table  I  lists 
the  M0RSE  module  routines  which  were  modified  to  permit  an  interface 


i 


I 


\ 

\  TABLE  I. 

•  \ 

M0HSE  Routines,  that  Depend  on  Which 
Geometry  Module  is  Used 

EUCLID  (3) 

I  .  G0MST  (2) 

INPUT  1 
MS0UR  (1) 

\  _ 

with  the!  CG  module.  The  numbers  after  the  routines  Relate  the  features 
mentioned  above  with  the  function  of  the  routines.  The  routines  without 
numbers  are  concerned  with  the  input  of  geometry  information.  Table  n 

TABLE  II. 


Geometry  Module  Routines 


05R  Spherical  Geometry 

Combinatorial  Geometry 

GE0M 

(2,3) 

GENI 

G0MFLP 

(5) 

GTVLIN 

J0MIN 

GG 

^  (2,  3) 

1  (2,3) 

'i  L00KZ 

(1) 

G1 

N0RML 

(4) 

ALBERT 

GTV0L 

\ 

PR 

G0MFLP 

J0MIN 

(5)  \ 

\ 

i 

L00KZ 

(1) 

N0RML 

GTV0L 

(4) 

2 


_ _ 


lists  the  geometry  module  routines  for  the  05R  spherical  and  the  com¬ 
binatorial  geometry  packages. 

Only  minor  changes  were  made  to  INPUT  and  MS0UR  and  no 
additional  documentation  is  provided.  INPUT  was  split  into  two  routines, 
INPUT1  and  INPUT2,  to  allow  a  more  efficient  overlay.  Numerous 
changes  were  required  for  subroutine  EUCLID  and  G0MST  and  new  de¬ 
scriptions  and  flow  charts  of  these  routines  aye  given.  A  complete 
description  of  the  combinatorial  geometry  package  is  included  and 
sample  problems  are  provided. 
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II. 


COMBINATORIAL  GEOMETRY 


2. 1  Body  Types 

Combinatorial  geometry  (CG)  describes  general  three  dimen¬ 
sional  material  configurations  by  considering  unions,  differences 
intersections  of  simple  bodies  such  as  spheres,  boxes,  cylinders,  etc. 

In  effect,  the  geometric  description  subdivides  the  problem  space  into 
♦ 

unique  zones.  Each  zone  is  the  result  of  combining  one  or  more  of 
the  following  geometric  bodies. 

1.  Rectangular  Parallelepiped  (RPP) 

2.  Box  (An  RPP  randomly  oriented  in  space) 

3.  Sphere 

4.  Right  Circular  Cylinder 

5.  Right  Elliptic  Cylinder 

6.  Truncated  Right  Angle  Cone 

7.  Ellipsoid 

8.  Right  Angle  Wedge 

9.  Arbitrary  Convex  Polyhedron  of  4,  5,  or  6  sides. 

Body  types  S-9  may  be  arbitrarily  oriented  with  respect  to 
the  x,  y,  z  coordinate  axes  used  to  determine  the  space.  Body  1,  a 
special  body  described  below,  must  have  sides  which  are  parallel  to 
the  coordinate  axes. 

The  basic  technique  for  the  description  of  the  geometry  con¬ 
sists  of  defining  the  location  and  shape  of  the  various  zones  in  terms  of 

* 

To  avoid  confusion  between  importance  regions  and  combinatorial 
geometry  regions,  we  depart  from  previous  combinatorial  geo¬ 
metry  descriptions  and  use  the  term  zone  to  indicate  a  combinatorial 
geometry  region  which  is  designated  by  the  variable  IR.  The  term 
region  is  reserved  for  an  importance  region.  Thus  the  zone  index 
is  IR. 


the  intersections  and  unions  of  the  geometric  bodies.  A  special  oper¬ 
ator  notation  involving  the  symbols  (+),  (-),  and  (0R)  is  used  to  describe 
the  intersections  and  unions.  These  symbols  are  used  by  the  program 
to  construct  information  relating  material  descriptions  to  the  body  de¬ 
finitions. 

If  a  body  appears  in  a  zone  description  with  a  (+)  operator,  it 
means  that  the  zone  being  described  is  wholly  contained  in  the  body.  If 
a  body  appears  in  a  zone  description  with  a  (-)  operator,  it  means  that 
the  zone  being  described  is  wholly  outside  the  body.  If  the  body  appears 
with  an  (0R)  operator,  it  means  that  the  zone  being  described  includes 
all  points  in  the  body.  In  some  instances,  a  zone  may  be  described  in 
terms  of  subzones  lumped  together  by  (0R)  statements.  When  (0R) 
operators  are  used  there  are  always  two  or  more  of  them,  and  they  re¬ 
fer  to  all  body  numbers  following  them,  either  (0  or  (-). 

Techniques  for  describing  a  particular  geometry  are  best  il¬ 
lustrated  by  examples.  Consider  an  object  composed  of  a  sphere  and  a 
cylinder  as  shown  in  Fig.  1.  To  describe  the  object,  one  takes  a  spher¬ 
ical  body  (2)  penetrated  by  a  cylindrical  body  (3)  (see  Fig.  1).  If  the 
materials  in  the  sphere  and  cylinder  are  the  same,  then  they  can  be 
considered  as  one  zone,  say  zone  I  (Fig.  lc).  The  description  of  zone  I 
would  be 

I  =  0R  +  2  0R  +  3 

This  means  that  a  point  is  in  zone  I  if  it  is  either  inside  body  2  or  inside 
body  3. 


If  different  materials  are  used  in  the  sphere  and  cylinder, 
then  the  sphere  with  a  cylindrical  hole  in  it  would  be  given  a  different 
zone  number  (say  J)  from  that  of  the  cylinder  (K). 

The  description  of  zone  J  would  be  (Fig.  Id): 

J  =  +  2  -  3 

This  means  that  points  in  zone  J  are  all  those  points  inside  body  2 
which  are  not  inside  body  3. 

The  description  of  zone  K  is  simply  (Fig.  le): 

K  =  +  3 


That  is,  all  points  in  zone  K  lie  inside  body  3. 

Combinations  of  more  than  two  bodies  and  similar  zone  de¬ 
scriptions  could  contain  a  long  string  of  (+),  (-),  and  (0R)  operators. 

It  is  important  however  to  remember  that  every  spatial  point  in  the 
geometry  must  be  located  in  one  and  only  one  zone. 

As  a  more  complicated  example  of  the  use  of  the  (0R)  operator, 
consider  the  system  shown  in  Fig.  2  consisting  of  the  shaded  zone  A  and 
the  unshaded  zone  B.  These  zones  can  be  described  by  the  two  B0X's, 
bodies  1  and  3,  and  the  RCC,  body  2.  The  zone  description  would  be 


and 


A  =  +  1  +  2 

B  =  0R  +3-1  0R  +3-2 


Notice  that  the  0R  operator  refers  to  all  following  body  numbers  until 
the  next  0R  operator  is  reached. 
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Fig.  2.  Use  of  0R  operators. 

The  geometry  must  be  specified  by  establishing  two  tables. 
The  first  table  describes  the  type  and  location  of  the  set  of  bodies  used 
in  the  geometrical  description.  The  second  table  identifies  the  physical 
zones  in  terms  of  these  bodies.  The  input  routine  processes  these 
tables  to  put  the  data  in  the  form  required  for  ray  tracing.  Because  the 
ray  tracing  routines  cannot  track  across  the  outermost  body,  all  of  the 
zones  must  be  within  a  surrounding  external  void  so  that  all  escaping 
particles  are  absorbed.  Also  no  point  may  be  in  more  than  one  zone. 

The  information  required  to  specify  each  type  of  body  is  as 

follows: 

a.  Rectangular  Parallelepiped  (RPP) 


Specify  the  minimum  and  maximum  values 
of  the  x,  y,  and  z  coordinates  which  bound 
the  parallelepiped. 


z 


Fig.  3.  Rectangular  Parallelepiped  (RPP). 
b)  Sphere  (SPH) 

Specify  the  vertex  V  at  the  center  and  the  scalar, 
R,  denoting  the  radius. 


c)  Right  Circular  Cylinder  (RCC) 

Specify  the  vertex  V  at  the  center  of  one  base, 
a  height  vector,  H ,  expressed  in  terms  of  its 
x,  y,  and  z  components,  and  a  scalar,  R, 
denoting  the  radius. 


„Ymax 
— Zmax 
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Fig.  5.  Right  Circular  Cylinder  (RCC). 

d)  Right  Elliptical  Cylinder  (REC) 

Specify  coordinates  of  the  center  of  the  base 
ellipse,  a  height  vector,  and  two  vectors  in 
the  plane  of  the  base  defining  the  major  and 
minor  axes.  Presently  this  body  is  not 
implemented. 


Fig.  6.  Right  Elliptical  Cylinder  (REC). 
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Truncated  Right  Angle  Cone  (TRC) 

Specify  a  vertex  V  at  the  center  of  the 
lower  base,  the  height  vector,  H,  expressed 
in  terms  of  its  x,  y,  z  components,  and  two 
scalars,  Rj  and  R2  ,  denoting  the  radii  of  the 
lower  and  upper  bases. 


Fig.  7.  Truncated  Right  Angle  Cone  (TRC). 
Ellipsoid  (ELL) 


Specify  two  vertices,  VjandVg,  denoting 
the  coordinates  of  the  foci  and  a  scalar,  R, 
denoting  the  length  of  the  major  axis. 


F 

1 


Specify  the  vertex  V  at  one  of  the  corners 
by  giving  its  (x,  y,  z)  coordinates.  Specify 
a  set  of  three  mutually  perpendicular 
vectors,  ai  ,  with  a^  and  ^  describing  the 
two  legs  of  the  right  triangle  of  the  wedge. 
That  is,  the  x,  y,  and  z  components  of  the 
height,  width,  and  length  vectors  are 
given. 

h 

v 

Fig.  9.  Right  Angle  Wedge  (WED). 

Box  (B0X) 

Specify  the  vertex  V  at  one  of  the  corners  by 
giving  its  (x,  y,  z)  coordinates.  Specify  a  set 
of  three  mutually  perpendicular  vectors,  a^ 
representing  the  height,  width,  and  length  of 
the  box,  respectively.  That  is,  the  x,  y, 
and  z  components  of  the  height,  width,  and 
length  vectors  are  given. 
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Vx,  Vy,  Vz 

Fig.  10.  Box  (B0X). 

i)  Arbitrary  Polyhedron  (ARB) 

Assign  an  index  (1  to  8)  to  each  vertex.  For 
each  vertex,  give  the  x,  y,  z  coordinates. 
Each  of  the  six  faces  are  then  described  by  a 
four -digit  number  giving  the  indices  of  the 
four  vertex  points  in  that  face.  For  each  face 
these  indices  must  be  entered  in  either  clock¬ 
wise  or  counterclockwise  order. 

Ys 


Fig.  11.  Arbitrary  Polyhedron  (ARB),' 
Detailed  input  instructions  are  given  in  Appendix  B. 
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2.2 


Subroutines 


Descriptions  of  the  combinatorial  geometry  routines  are  given 
with  logical  flow  charts  depicting  their  functions.  Two  M0RSE  routines, 
G0MST  and  L00KZ  required  significant  modifications  for  use  with  the 
CG  package  and  descriptions  of  these  routines  are  also  provided. 

2.  2. 1  Subroutine  G1(S,  MA,  FPD,  L0CREG,  NUMB0D,  IR0R, 

ml,  IR2) 

G1  is  the  control  routine  for  the  combinatorial  geometry.  On 
one  call,  it  calculates  the  distance  travelled  in  the  present  zone,  and 
the  number  IR  of  the  next  zone  to  be  enterec*.  Essentially  GG  is  called 
for  each  body  adjacent  to  the  present  zone,  calculating  RIN  and  R0UT, 
the  distances  to  entry  and  exit  of  the  body  along  the  trajectory.  The  next 
zone  to  be  entered  is  determined  by  again  calling  GG  to  calculate  RIN 
and  R0UT  for  each  body  adjacent  to  the  next  possible  zone.  These  next 
possible  zones  are  determined  by  examining  a  list  cf  all  the  previous 
zones  entered  on  crossing  this  body.  RIN  and  R0UT  are  checked  against 
the  input  zone  descriptions  to  determine  the  correct  zone.  If  it  is  not 
found  in  the  list  of  previous  zones,  all  other  zones  are  examined  in  a 
similar  fashion,  and  when  the  correct  zone  is  found,  it  is  added  to  the 
list  of  previous  zones  for  that  body.  If  the  new  zone  is  different  from 
the  old,  G1  returns;  otherwise  G1  continues  tracking  until  a  different 
zone  is  encountered.  One  change  added  to  the  M0RSE  version  of  G1  is 
that  if  the  distance  to  the  next  boundary  is  greater  than  the  distance  to 
scattering,  G1  returns  without  determining  the  next  zone  past  the  bound¬ 
ary,  setting  the  flag  MARKG  in  common  0RGI. 

Called  from:  G0MST,  EUCLID,  MESH 

Subroutines  called:  GG 

Commons  required:  PAREM,  G0ML0C,  DBG,  0RGI 
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Variables 


\ 

required:  ! 

XB(3) 

WB(3) 

IR 

DIST 

DISTO 

NASC 

KL00P 

PINF 

MA,  FPD, 
IR1,  IR2 


■  i 

starting  coordinates  of  present  trajectory, 
direction  cosines  of  trajectory,  ^ 

present  zone,  ® 

present  distance  from  XB(3), 

distance  frbm  XB(3)  to  next  scattering  point, 

I  !  \ 

less  ihan  zero  if  this  is  a  new  trajectory, 

trajectory  index, 
machine  infinity, 

L0CREG,  NUMB0D,  IR0R,  ' 

locations  in  blank  common  used  for  variable 
dimensioning. 


Variables  changed: 


KL00P 

NASC  ' 
LSURF 

\  ,  . 

I  DIST  ' 

MARK 

S 

IRPRIM 
MA(INEXT)  - 

MA(INEX)  - 


trajectory  index  incremented  if  this  is  a 
new  trajectory, 

next  fyody  intersected  by  trajectory, 

surface  of  body  NASC  crossed  at  next  inter¬ 
section  (negative  if  leaving  and  positive  if 
entering  NASC),  ' 

distance  from  XB(3)  to  next  intersection  or 
collision  site,  .  » 

\ 

set  to  1  if  distance  to  collision!  (DISTO)  is 
less  than  distance  to  next  intersection  (other¬ 
wise  0),  1 

distance  travelled  on  this  call  to  Gl, 

zone  to  be  entered  on  boundary  crossing, 

new  zone  added  to  the  list  of  next  possible 
zones  for  body  NB0, 

location  in  MA  of  next  item  in  the  next  possible 
zone  list  for  body  NB0  (these  lists  leap-frog 
through  the  end  of  the  MA  array). 


i 
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Significant  internal  variables: 


NB0 

\ 

R0UT 

RIN 

LRI 
LR0 
!  IRP 

Limitations: 

I0UT 


absolute  value  is  body  being  considered 
white  a  negative  or  positive  sign  indicate^ 
that  zone  IR  or  IRP  is  outside  or  inside  the 
body  respectively, 

distance  to  exit  of  body  NB0  calculated  by 
GG, 

distance  to  entry  of  body  NB0  calculated  by 
GG,  , 

surface  of  body  NB0  entered  by  trajectory, 

surface  of  body  NB0  trajectory  exits, 

zone  being  considered  as  next  zone. 


output  logical  unit  is  set  to  6  in  data  state¬ 
ment. 


1 
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Subroutine  G1  (S,  MA,  FPD,  L0CREG,  NUMB0D,  IH0R,  IR1,  IR2) 


Get  next  possible  zone 


Is  next  possible  zone 
IRP  zero? 


Loop  over  all  bodies  NB0  defining  zone 
IRP  to  determine  if  ERP  is  next  zone. 


Call  GG 


Yes— 


Is  NB0  negative?  ) - ►No 


'is  R0UT  <  0  or\_vjr 
R0UT  <  DIST  r 
lor  RIN  > DIST ?/ 


Not  this 


Is  RIN<  DIST) 
and 

DIST  <  R0UT?J 


Found  a  zone  j 


Set  IRPRIM 


If  IRPRIM  =  IR - -  No  RETURN 


Move  point  to  next  boundary 
and  begin  again. 


:<s 


2. 2. 2  Subroutine  L00KZ  (X,  Y,  Z,  MA,  FPD,  L0CREG, 

NUMB0D,  IR0R,  NS0RT 

The  purpose  of  this  routine  is  to  return  the  combinatorial  geo¬ 
metry  zone  of  point  (X,  Y,  Z)  so  that  tracking  can  be  initialized.  The 
coding  has  been  borrowed  from  the  second  half  of  subroutine  G1  and 
adapted  to  determine  the  zone  of  a  source  particle.  For  efficiency 
L00KZ  builds  a  list  of  possible  source  zones  to  search  on  future  calls. 

If  the  region  is  not  found  on  this  list,  all  other  zones  are  examined  and 
upon  determining  the  new  source  zone,  it  too  is  added  to  the  list.  Notice 
that  the  starting  direction  cosines  (.8,  .6,  0. 0)  are  assumed  in  L00KZ, 
but  may  be  changed  elsewhere. 


Routines  called:  GG 

Commons  required:  0RGI,  PAREM,  G0ML0C,  and  DBG. 
Variables  required: 


X,  Y,  Z  -  coordinates  for  which  a  zone  is  desired, 

MA,  FPD,  L0CREG,  NUMB0D,  IR0R,  NS0R 

locations  in  blank  common  used  for  variable 
dimensioning. 


Varaibles  changed: 

KL00P  -  trajectory  index  is  incremented, 

NMED  -  common  0RGI  variable  set  to  correct  zone 
number, 


NS0R(INEXT)  -  new  zone  added  to  list  of  possible  source 
zones. 


Significant  internal  variables: 

WB(3)  -  set  to  .  8,  .  6,  0. 0  so  that  L00KZ  need  not 

be  called  with  a  direction 


♦ 

CG  version. 
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Subroutine  L00KZ  (X,  Y,  Z,  MA,  FPD,  L0CREG,  NUMB0D,  IR0R, 

NS0R) 


2. 2. 3  Subroutine  GG  (L0CAT,  MA,  FPD) 

GG  is  the  workhorse  of  the  combinatorial  geometry,  comput¬ 
ing  distances  to  intersections  for  all  body  types.  It  is  called  from  G1 
or  L00KZ  to  compute  distance  to  entry  and  distance  to  exit  to  a  body 
whose  location  in  the  MA  array  is  flagged  by  the  argument  L0CAT.  Each 
time  a  distance  to  entry,  RIN,  and  a  distance  to  exit,  R0UT,  are  calcu¬ 
lated  for  a  body  they  are  stored  in  the  MA  array  together  with  LRI  and 
LR0,  the  indices  of  the  entry  and  exit  surfaces  of  the  body.  Also  stored 
at  this  time  is  KL00P,  the  particles  trajectory  index.  On  a  subsequent 
call  to  GG  for  that  body,  KL00P  is  checked  against  the  earlier  value, 
now  L00P.  If  they  are  the  same,  the  old  values  of  RIN,  R0UT,  LRI, 
and  LR0  are  retrieved  so  that  GG  can  return  immediately. 

If  it  is  necessary  to  compute  a  new  trajectory  a  different  area 
of  coding  is  entered  for  each  body  type  to  calculate  RIN,  R0UT,  LRI 
and  LR0. 

Called  from:  G1  and  L00KZ 

Subroutines  called:  none 
Commons  required:  PAREM,  GE0M2 
Variables  required: 

LOCAT  -  starting  locati  n  in  the  MA  array  of  integer 
data  for  the  appropriate  body, 

KL00P  -  trajectory  index, 

PINF  -  machine  infinity  (stored  in  RIN  and  R0UT 

when  the  trajectory  misses  the  body), 

MA,  FPD  -  locations  in  blank  common  required  for 
variable  dimensioning. 
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Variables  changed: 


RIN 

distance  to  entrance, 

R0UT 

distance  to  exit, 

LRI 

-  surface  of  entrance 

LR0 

surface  of  exit. 

Significant  internal  variables: 

L  -  same  as  L0CAT,  starting  location  in  MA 

array  of  body  data  for  the  appropriate  body, 

K  -  starting  location  in  FPD  array  of  floating 

point  data  for  the  appropriate  body. 
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Subroutine  GG  (L0CAT,  MA,  FPD) 


START 


Retrieve  K,  L00P,  and 
IT  YPE  from  MA. 

_ _  I  Set  LRI,  LR0,  RIN 

Has  this  trajectory  previously  T  and  R0UT  from  MA 


been  calculated  for  this  body?. 


|and  FPD  arrays. 


L  Is  IT  YPE  valid? 


-no*  Write  error  (RETURN) 
message, _  V  J 


Store  KL00P  in  MA  array! 
_ for  this  body. _ | 


(STOP) 


Transfer  on  IT  YPE  to  compute  RIN,  R0UT,  LRI,  and  LR0 


2.2.4 

Although  subroutine  J0MIN  is  used  to  read  in  all  of  the  com¬ 
binatorial  geometry  data  except  the  region  volumes,  that  data  is  put  into 
the  proper  storage  location  in  blank  common  by  GENI.  This  is  accom¬ 
plished  by  having  J0MIN  write  the  body  and  zone  data  on  a  mass  storage 
unit,  so  that  GENI  can  retrieve  that  data.  GENI  assigns  the  data  to 
blank  common  in  the  area  set  aside  by  J0MEN.  For  an  arbitrary  poly¬ 
hedron  body,  GENI  calls  ALBERT  to  handle  a  portion  of  the  data.  In 
addition,  GENI  outputs  the  input  information  concerning  the  geometry 
data.  GENI  computes  certain  geometry  dependent  data  needed  during 
the  random  walk. 

Called  from:  J0MIN 


Subroutines  called: 

ALBERT 

Commons  required: 

Blank,  G0ML0C,  PAREM 

Variables  required: 

11 

input  unit, 

I0UT 

output  unit, 

IN 

bulk  storage  logical  unit. 

Variables  input:  Geometry  data  for  zones  and  bodies  is  retrieved 

from  mass  storage  unit  IN. 

Variables  changed:  MA(I),  FPD(I),  L0CREG(I),  IR0R(I),  MRCZ(I), 

MMCZ(I),  KR2(I)--see  Fig.  13,  layout  of  com¬ 
binatorial  geometry  storage  in  blank  common. 


Subroutine  GENI  (MA,  FPD,  L0CREG,  NUMB0D,  IR0R, 
MRIZ,  MRCZ,  MMIZ,  MMCfc,  KR1,  KR2,  il,  I0UT,  IN) 


Subroutine  GENI  (MA,  PPD,  L0CREG,  NUMB0D,  IR0R,  MRIZ,  MRCZ 
MMIZ,  MMCZ,  KR1,  KR2,  IN) 


2.  2.  5  Subroutine  ALBERT  (F,  IERR) 

ALBERT  is  called  by  subroutine  GENI  to  process  the  arbi¬ 
trary  polyhedron  (ARB)  body  data  before  storage  in  the  FPD  array.  The 
ARB  body  data  as  read  by  GENI  consists  of  the  coordinates  of  all  corners 
(eight  for  a  six-sided  figure),  followed  by  a  packed  decimal  number  for 
each  side  indicating  which  corner  points  make  up  that  side.  ALBERT 
processes  this  data,  returning  a  unit  normal  vector  and  a  minimum 
distance  to  origin  for  each  plane  containing  a  side  of  the  ARB.  An  ARB 
can  have  up  to  six  sides.  These  unit  vectors  and  distances  then  replace 
the  original  ARB  data  in  the  FPD  array.  The  number  of  sides  and  a 
distance  that  is  characteristic  of  the  ARB’s  minimum  dimension  are 
also  stored  in  the  FPD  array.  This  minimum  distance  is  later  ust ...  for 
round-off  tests. 

Called  from:  GENI 

Subroutines  called:  none 

Commons  required:  none 

Variables  required: 

F  -  FPD  array  read  from  GENI 

Variables  changed: 

F  -  FPD  array  changed  to  contain  unit  normals 

and  distances  to  origin  for  each  plane.  Num¬ 
ber  of  sides  and  a  minimum  distance  are 
also  stored. 

Significant  internal  variables: 

X(8, 3)  -  coordinates  of  cornei  points, 

IS(6, 4)  -  indicates  corners  contained  in  each  of  the 

possible  six  surfaces, 

V(3, 4)  -  edge  vectors  for  a  given  surface  (an  inner 

cross  product  then  gives  the  normal  vector). 
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Subroutine  ALBERT  (F,  IERR) 


(START) 


[store  coordinates  of  corner  points  in  X(8, 3)  and  zero  F  array, 
(Count  the  sides  and  store  corner  indices  in  IX(6, 4).| 


Calculate  minimum  distance  between  points.] 


» _  ... 

["Begin  a  calculation  for  each  side. 

|  ' 

[Compute  two  adjacent  edge  vectors. 

_ LL _ 

Compute  normal  vectors  as  the  cross  products 
_ of  two  adjacent  edge  vectors. _ 


[Compute  niinimum  distance  to  orgin  as  dot  product  of 
l  normal  vector  and  coordinates  of  the  corner. 
- ] - 

, - * - ; - 

Count  the  number  of  points  on  either  side  of  the  plane. 


Are  all  points  on  the  far  side  of  the  plane?) - Yes 


I  Charge  sign  of  the  normalizing  factor.  | 

Normalize  normal  vectors  to  unity.  Also 
normalize  the  distance  by  the  seme  factor. 


_L 


1  Is  this  the  last  side  ?  } -  no - 

yes' 

Store  normal  vectors,  distances,  minimum  side  lengtn 
and  numocr  of  sides  in  FPD  array. 

.^RETURN\ 
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Subroutine  PR  (K) 


32 


2.2.7  Subroutine  N0RML  (MA,  FPD,  L0CREG,  NUMB0D)* 

The  purpose  of  subroutine  N0RML  is  to  return  a  unit  vector 
to  a  combinatorial  geometry  body  NASC  at  point  XB  +  D  1  *  WB  which 
must  be  on  the  surface  LSURF  of  body  NASC.  This  unit  vector  is  useful 
either  for  albedo  scattering  or  boundary  crossing  flux  estimates.  The 
sign  of  the  unit  vector  is  chosen  so  that  wB  •  ri  is  negative  meaning  that 
the  unit  vector  will  point  against  the  particle  direction. 

Called  from:  ALBD0 

Commons  required:  N0RMAL,  G0ML0C,  PAREM 

Variables  required: 

NASC  -  body  number  particle  is  on, 

LSURF  -  surface  of  body  NASCr 

XB(3) 

WB(3)  -  trajectory  parameters, 

DIST 

MA,  FPD,  L0CREG,  NUMB0D 

locations  in  blank  common  used  for  variable 
dimensioning. 

Variables  changed: 

UN(3)  -  direction  cosines  of  unit  vector  stored  in 

common  N0RMAL. 

Significant  internal  variables: 

XP(3)  -  set  to  XB  +  DIST  *  ifiTJ, 

X(3)  -  usually  used  to  relate  the  intersection  point 

to  a  body  centered  coordinate  system, 

H(3)  -  usually  used  as  a  body  orientation, 

L  -  body  location  in  MA  array, 

K  -  body  location  in  FPD  array. 

It 

CG  version. 
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2.2.8 


Subroutine  r~VLIN  (FPD,  MRIZ,  VN0R,  IV0PT,  NIR, 

niz,  ii,  10;  - - - - — 

This  routine  is  called  by  J0MIN  to  read  in  or  to  calculate  the 
volume  of  each  region  in  the  geometry.  The  four  options  available  are 
(1)  to  set  each  volume  to  1. ,  (2)  to  calculate  the  volumes  for  concentric 
spheres,  (3)  to  calculate  the  volumes  for  slabs  (not  coded  at  present), 
or  (4)  to  read  in  the  volumes  for  each  region  from  cards.  The  volumes 
are  stored  in  blank  common  and  are  only  used  by  the  track  length  and 
collision  density  estimators. 

Called  from:  J0MIN 

Subroutines  called:  none 
Commons  required:  none 
Variables  required: 

FPD  -  array  containing  zone  description  data, 

MRIZ  -  array  which  relates  M0RSE  region  to  input 

zones, 

VN0R  -  array  containing  the  M0RSE  volumes, 

IV0PT  -  options  for  determining  volumes, 

NIR  -  number  of  regions, 

NIZ  -  number  of  zones, 

II,  10  -  input  and  output  logical  units. 

Variables  changed: 

Variables  in  blank  common  starting  at  KV0L. 


2.2.9  Subroutine  J0MIN  (NADD,  II,  IQ)* 

The  input  of  the  geometry  data  is  controlled  by  the  J0MIN 
subroutine,  which  performs  the  following  tasks: 

•  Reads  all  geometry  input  data  except  the  region  volumes. 

•  Writes  the  body  and  zone  data  on  a  mass  storage  unit 
(I0UT-16). 

•  Determines  the  length  of  all  geometry  arrays. 

•  Calculates  the  beginning  location  in  blank  common  of 
geometry  arrays. 

•  Initializes  geometry  arrays. 

•  Calls  the  GTVLIN  subroutine  which  returns  region 
volumes. 

Since  combinatorial  geometry  input  data  is  dynamically  allocated  to  con¬ 
serve  storage  area,  it  is  stored  temporarily  on  a  mass  storage  device. 
This  allows  the  core  storage  requirements  to  be  determined.  Hence, 
much  of  the  coding  in  J0MIN  is  similar  to  GENI,  which  reads  the  data 
on  the  mass  storage  device  and  puts  it  into  blank  common. 

Called  from:  INPlfol 

Subroutines  called:  GENI,  GTVLIN 
Commons  required:  B^ank,  G0ML0C,  PAREM 
Variables  required:  ^ 

4 

All  variable^'  in  G0ML0C,  II,  10 
Variables  input:  : 

IV0PT,  IDBp,  NAZ,  MRIZ(I),  MMIZ(I) 

A 

Variables  changed:  A 

A 

All  variable^  in  G0ML0C,  IV0PT,  IDBG,  MRIZ(I),  MMIZ(I). 

A 

*  A 

CG  version. 
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Important  internal  variables: 


I0UT 

NAZT 


mass  storage  unit, 

total  number  of  adjacent  zones  summed  over 
all  zones. 


38 


■tfiHrtl  MrYiiil'initlri 


Subroutine  J0MIN  (NADD,  11,  10) 
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2 . 2. 1°  Subroutine  GTV0L  (MXREG,  GN0R) 

This  routine  is  called  by  ENDRUN  to  calculate  the  reciprocal 
of  the  volume  of  each  region.  It  is  geometry  independent,  since  it  uses 
information  determined  in  GTVLIN. 

Called  from:  ENDRUN 

Subroutines  called:  none 
Commons  required:  G0ML0C,  blank 
Variables  required: 

KV0L  -  index  from  common  G0ML0C, 

MXREG  -  number  of  regions  in  the  geometry. 

Variables  changed: 

GN0R  -  reciprocal  of  the  volume  io r  e  ion. 


2.2.11  Subroutine  G0MST  (TSIG,  MARK)* 

Any  boundary  crossing  between  the  present  and  next  collision 
sites  are  determined  by  calling  the  combinatorial  geometry  routine  Gl. 
Before  the  Gl  call,  combinatorial  geometry  variables  in  common  PAREM 
are  initialized,  and  after  the  call  NUTR0N  variables  are  updated. 

WATE*S  is  summed  in  blank  common  as  a  track  length  flux  estimator. 
MARK  is  set  to  -1  for  an  external  void  and  -2  for  an  internal  void. 


Called  from: 
Subroutines  called: 
Commons  required: 
Variables  required: 


NXTC0L 

Gl 

NUTR0N,  AP0LL0,  PAREM,  0RGI,  G0ML0C 


X,  Y,  Z,  U,  V,  W,  NMED  see  common  NUTR0N  for  definitions. 


BLZNT 

XB(3) 

WB(3) 

MARK 

ETATH 

DIST 

Variables  changed: 

X,  Y,  Z 

DISTO 

ETAUSD 

ETATH 


value  of  IR  from  last  track  or  from  L00KZ, 

starting  coordinates  of  present  trajectory, 

trajectory  direction  cosines, 

flag  to  indicate  type  of  trajectory, 

distance  to  be  travelled  in  cm  if  the  flight 
remains  in  the  same  media, 

present  distance  from  XB(3). 


endpoints  of  flight, 

distance  from  XB(3)  to  next  collision  site, 

mean  free  paths  travelled  on  this  call  to 
G0MST, 

cm  travelled  on  this  call  to  G0MST, 


* 

CG  version. 
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MARK  -  flag  indicating  type  of  termination  of  flight, 

0  -  normal  boundary  crossing, 

1  -  flight  within  one  medium , 

-1  -  particle  escaped, 

-2  -  particle  entered  an  interior  void, 

NMED  and 

NREG  -  medium  and  region  of  end  point, 

combinatorial  geometry  region  of  end  point. 


BLZNT 


2. 2. 12  Subroutine  EUCLID  (MRK,  XI,  Yl,  Zl,  X2,  Y2,  Z2,  P1P2, 
IG)  ARG,  NT,  MEDIUMTIBLZ,  NREGN)* 

'  \  \ 

Subroutine  EUCLID  is  used  to  determine  the  number  of  mean 
free  paths  between  two  points.  However,  if  the  argument  NT  is  non¬ 
zero  only  the  mean  free  'paths  to  the  first  media  boundary  will  be  calcu¬ 
lated.  The  routine  works  much  like  NXTC0L  and  G0MST.  After  ini- 

tializihg  the  combinatorial  geometry  variables,  the  particle  is  tracked 

'  \  \ 

from  point  one  to  point  two.  Total  cross  sections  are  determined  by 
calling  NSIGTA  and  distances  to  next  interfaces  are  calculated  by  call¬ 
ing  Gl.  On  an  initial  call  to  EUCLID,  MRK  should  be  1.  If  NT  is  non¬ 
zero,  MRK  should  be  0  on  successive  calls  for  the  same  trajectory,  so 
that  new  trajectory  parameters  are  not  initialized.  ,  This  will  be  handled 
automatically  if  the  calling  routine  does  not  change  MRK. 


Called  from:  I 

GETETA,  RELC0L- 

Subroutines  called: 

Gl, 

NSIGTA 

Commons  required: 

PAREM,  0RGI 

Variables  required: 

\ 

\ 

MRK 

- 

set  to  1  for  initial  call  of  a  trajectory; 

SI,  Yl,  Zl 

- 

coordinates  of  starting  point, 

X2,  Y2,  Z2 

- 

coordinates  of  end  point, 

P1P2 

- 

distance  between  starting  and  end  points. 

IG  ' 

- 

energy  group  index, 

NT 

- 

0  for  total  mean  free  paths,  ^  0  for  mean 
free  paths  between  intersection  points, 

MEDIUM 

- 

NMED  of  point  one, 

NREGN 

- 

NREG  of  point  one,  , 

i 

IBLZN 

- 

IR  of  point  one. 

CG  version. 
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Variables  changed: 


MRK 


XI,  Y1.Z1 

ARG 

NT 

MEDIUM 


1  for  a  flight  reaching  the  end  point, 

0  for  a  flight  crossing  a  medium  boundary 
(NT  4  0  only), 

-1  for  a  flight  escaping  the  system, 

-2  for  a  flight  encountering  an  internal  void 
(NT  4  0  only), 

returns  boundary  intersection  point  if 
NT  4  0, 

negative  of  number  of  mean  free  paths, 

if  NT  4  0  on  input,  will  return  as  -1  if  an 
escape  occurs, 

medium  number  of  end  point. 


Significant  internal  variables: 

distance  remaining  to  point  two  in  cm, 
distance  travelled  on  last  call  of  Gl. 


ETA 

ETAUSD 


t» 


Subroutine  EUCLID 


RETURN 


2.  2. 13  Subroutine  ENDRUN 


This  subroutine  is  called  at  the  end  of  a  run.  The  weight  of 
the  particles  at  collision  and  the  track  length  of  the  particles  has  been 
summed  over  all  collisions  within  the  region  for  all  particles  in  the  run. 
Thus,  a  normalization  by  the  number  of  particles  in  the  run  is  required. 

In  addition,  the  volume  of  the  region  and  the  total  cross  section  for  each 
region  and  energy  group  is  used  to  determine  the  average  fluence  over 
the  regions  for  the  collision  density  estimate.  For  the  track  length  esti¬ 
mate  the  same  coding  is  used  with  the  total  cross  section  set  to  1.  Note 
that  a  region  cannot  contain  more  than  one  medium.  Since  storage  al¬ 
locations  are  flexibly  dimensioned,  NEX  in  the  SAMB0  input  must  provide 
for  MXREG  +  2  arrays  in  Blank  Common  to  be  used  by  ENDRUN.  A 
check  is  made  and  if  this  condition  is  not  met,  no  output  from  ENDRUN 
is  obtained. 

There  are  many  manipulations  in  ENDRUN  with  the  energy 
indexes.  These  manipulations  are  required  in  order  to  treat  the  six 
major  options  in  M0RSE.  Appendix  A  gives  an  example  of  the  six  prob¬ 
lems. 

Called  from:  NRUN 

Subroutines  called:  GTV0L,  NSIGTA 
Commons  required:  blank,  PDET,  USER 
Variables  required: 

NMED  (REGI0N)  -  must  be  set  in  data  statement, 

GN0R  -  reciprocal  of  the  volume  for  each 

region, 

Many  variables  in  blank  common. 


Significant  internal  variables: 

FNT  -  total  number  of  source  particles  in  the  run, 

MX REG  -  the  number  of  geometry  regions, 

LIM1  -  the  group  number  for  the  last  primary 

particle  group, 

EB0T1  -  lowest  primary  particle  group  energy, 

EB0T2  -  lowest  secondary  particle  group  energy, 

NGP3  -  the  group  number  for  last  secondary  group, 

NQP1  -  the  group  number  for  last  primary  group 

(negative  if  gamma-ray  only  problem), 

IFLG  -  flag  to  indicate  which  fluence  estimate  is 

currently  being  output.  0 -collision  density, 
1 -track  length. 


Subroutine  ENDRUN 
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O' 


I  <  NMGP 


Arc  there  more  primary 
_ particle  groups? 


Calculate  Etn(E)  for  each  secondary  group. 


III. 


MODIFICATIONS  FOR  UNIVAC-1108 


In  converting  the  M0RSE-SAMB0-M0RSEC  code  from  the 
IBM  360  to  the  UNIVAC-1108,  there  were  a  considerable  number  of 
changes  that  had  to  be  made.  These  changes  were  connected  mostly 
with  the  local  library  functions  of  the  IBM-360  and  an  attempt  was  made 
to  not  add  local  UNIVAC-1108  library  functions.  The  only  reduction  in 
the  code's  capability  is  in  the  diagnostic  module  where  the  ability  to 
look  at  a  bit  string  and  decide  whether  it  was  a  floating  point  or  integer 
number,  convert  it  to  hollerith  and  output  it,  does  not  exist  in  this  version. 
However,  the  diagnostic  module  is  mostly  machine  independent  and  still 
has  the  capability  of  writing  out  values  from  common  and  parts  of  blank 
common.  Subroutines  BNKHLP  and  HELPER  are  dummy  routines.  A 
diagnostic  routine  for  the  combinatorial  geometry  package  has  been  added. 

Routines  which  are  still  machine  dependent  (involving  coding 
other  than  word  lengths,  data  statements  and  equivalences)  are  DATE 
and  TIMER.  Subroutine  DATE  calls  ERTRAN  and  uses  DEC0DE.  Sub¬ 
routine  TIMER  c.  ils  CPUTIM  that  assumes  that  the  time  is  returned  in 
units  of  microseconds. 

Subroutine  READSG  is  still  machine  dependent  in  its  need  to 
read  each  card  twice  and  subroutine  ERR0R  has  been  added— its  only 
function  is  to  call  EXIT.  Subroutine  INPUT  now  calls  INPUT1  and 
INPUT2  with  this  change,  permitting  a  better  use  of  overlay.  The 
first  six  variables  in  NXTRA  from  common  AP0LL0  are  used. 

The  layout  of  blank  common  is  given  in  Figs.  12-18  for  the 
different  modules  of  M0RSE.  Definitions  of  variables  in  blank  common 
are  given  in  Table  m  which  defines  the  indexing  scheme.  Variables  in 
the  labelled  commons  are  defined  in  Tables  IV  -  XVI. 

Preceding  page  blank 


55 


A  description  of  the  energy  indexing  scheme  in  M0RSE- 
SAMB0-M0RSEC  is  given  in  Appendix  A. 

3. 1  Built-In  Fluence  Estimators 

There  are  two  fluence  estimators  that  are  "built-in”  the 
UNIVAC  1108  combinatorial  geometry  version  of  M0RSE.  The  addition 
of  the  collision  density  estimator  to  M0RSE  has  been  discussed  pre- 
viouslyv  and  the  track  length  per  unit  volume  estimator  has  been  added 
to  take  advantage  of  a  feature  of  the  combinatorial  geometry.  The  re¬ 
sulting  fluence  is  averaged  over  geometry  regions  and  if  no  regions  are 
used  only  the  spectrum  averaged  over  the  whole  system  is  available. 


3. 2  Collision  Density  Fluence  Estimator 


An  estimate  of  the  fluence  averaged  over  a  region  specified 
by  the  geometry  input  may  be  obtained  from  the  weight  that  has  been 
summed  over  all  collisions  in  the  region  by 


o 


V<Ei> 


Vi 

E 

n=l 


w 

n 


nstrt*nits*volv*£  Ty(E.) 


where 


fluence  in  region  V  at  energy  E., 


NITS*NSTRT 


number  of  particles  of  energy  E.  scatter¬ 
ing  in  region  V, 

weight  of  the  n  n  particle  of  energy  E. 
scattering  in  region  V, 

total  number  of  particles, 


VOLy  =  volume  region  V, 

W  ■ 
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tota1  cross  section  for  groups  E.  in  region 
V.  1 


3.3 


Track  Length  per  Unit  Volume  Estimator 


Another  estimate  of  the  fluence  averaged  over  a  region  may 
be  obtained  from  the  track  length  of  the  neutron  or  photon  in  that  region. 
The  track  length  calculated  by  the  combinatorial  geometry  package  is 
stored  in  a  manner  similar  to  the  collision  weights.  The  fluence  esti¬ 
mate  may  be  obtained  by 


where 


VEi> 


NVi 

T  L  w 
n=l  n  n 
NgTRTWT^V'OI, 


y 


NVi  =  number  of  trajectories  of  energy  E.  in  region  V, 

til 

wn  =  weight  of  particle  on  n  trajectory,  and 

til 

L  =  track  length  of  n  trajectory  in  region  V  at  energy 

n  E.  . 

1 


Subroutine  ENDRUN  is  called  to  process  the  weights  and 
track  lengths  to  obtain  estimates  of  the  fluence.  In  an  adjoint  problem, 
the  importance  is  calculated.  The  fluence  estimate  has  units  of  particles/ 
eV/cm  /source  and  in  an  adjoint  problem  the  units  are  importance/cm  / 
source.  The  total  fluence  is  integrated  over  energy  from  the  primary 
and  secondary  energy  groups  to  obtain  an  integral  value  for  each  region. 
The  only  input  required  for  this  routine  is  a  data  statement  giving  the 
medium  number  for  each  region;  there  is  a  limitation  of  only  one  medium 
^gion.  In  addition,  one  of  the  several  options  for  inputing  or  calcu¬ 
late.  .he  volume  of  each  region  must  be  utilized. 

The  major  disadvantage  with  using  these  estimators  is  that 
no  estimate  of  statistics  is  available.  The  use  of  this  estimator  is  in¬ 
expensive,  however  (output  routine  is  called  only  once),  and  is  meant 
to  be  used  as  a  supplemental  estimator. 
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Name 

Length 

KNETt 

VEL 

!  4*NMTG 

FS 

L0CWTS 

RFR  _ 

Current  Wei 
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>h{ 

|  3+MGPREG 

Path 

|  MGPRFG 

SPLIT  and  R.R. 
Counters 

J  8*MGPREG 

L0CFWL _ .. 

Initial  Weight 
Standards 

J  3*MGPREG 

Current  FWL0 

J  2*MXREG 

Initial  FWLD 

L.0CEPR 

GW  1.0 

|  NMTG*MXREG 

LOCNSC 

EPRB 

}  NMTG*MXREG 

Scattering 

Counters 

}  10+NMTG+MXREG 

L0CrSN  _ _ 

NSCA 

|  MEDIA 

FISH 

|  NMTG*MEDIA 

FSE 

j  NMTG*  MEDIA 

NGE0M 

GMGN 

|  NMTG+  MEDIA 

nglast _ * 

Geometry  Data 

j  See  Figs.  13-15  for 
details. 

NS1GL 

Permanent  Cross 
Sections 

|  See  Fig.  16  for  details. 

NMXsrcr:  ...  _  .  r 

Temporary 

Cross 

Sections 

Particle 

Bank 

1  1  INMOST  see  Fig.  17 
|  for  details. 

NLAST 

Fission 

Bank 

\  7  "N  MO  ST  IF  MFISTP>0 
j  see  Fig.  18  for  details. 

User  Arcv1 

1  NLEFT  son  Fig.  1 9  for 
j  details. 

Fig.  12.  General  layout  of  blank  common 
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\ 
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KMCZ 

KKR1 

KKR2 

lO-ISR 


MA  i 
|  Integer  array  j 


FPD 

Floating  point  array 


L0CREG 

Indices  to  correlate  MA  array 
data  wi’n  code  zone  data  ■ 

NU1^B0D 

Number  of  bodies  for  each 
code  zone 

IR0R 

Indices  to  correlate  input  zone 
to  code  zone  i 

- T - 1 - \ 

MRIZ  j 

Indices  to  correlate  M0RSE 
region  to  input  zone 

MRCZ 

Indices  to  correlate  M0RSE 
region  to  code  zone 

NM3Z 

Indices  to  correlate  M0RSE 
media  to  input  zone  I 
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zone  to  ‘  nput  zones 

KR2 

Indices  to  correlate  last  code 
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KV0L 
NGIiAST  % 


NS0R 

Indr-cs  of  code  zones  in  which 
source  particles  have  been  found 
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\  Volume  of  each  M0RSF  region 
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NUMR 

1 
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IRTRV 
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IRTRV 
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Fig.  13.  Layout  of  combinatorial  geometry  data  in  bla  *xiv  common. 
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Fig.  14.  Detailed  layout  of  the  FPD  array  in  blank  common. 
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Figure  15. 
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Size 


1RSG  ■-  NLAST 
=  NGI.AST+1 

INGS 

INSG 

IDEL 


List  of  Mixing  Table 


Index  to  Egg  (Primary) 

Index  to  Egg  (Secondary) 

List  of  Element  I.D.  Numbers 


INNN 

IGGG 

ISTART 

JSCC0G 

INAB0G 

IGAB0G 

IFP0RG 

INUS 

INGN 

INGNP 

IFNGT> 

IFSP0G 

IDSG0G 

IPRBNG 

IPRBGG 

ISCANG 

ISCAGG 

ISP0RG  +  ISTART 
If  ISTAT  >  0  ISP0RT 


NLAST  or  NNIC 
NNIC 


Number  of  Downscatters  for  each 
primary  group 

Number  of  Downscatters  for  each 
secondary  group _ 


as/LT _  _ _ _ ' 

V£  T _ 

vLf/£T 

V~g _ 

En/Et  not  used  at  present 


E  _  not  used  at  present 

£N-r _ 

EN, 

g-g _ 


E  y 


fLriL 


PN, 
g  ~B 


p?. 


B  -8 


AN, 

V~B 


6  -6 

Repeat  for  next  medium 


l'i  Coefficient  Primary 
Pj  Coefficient  Secondary 
Repeal  for  Coefficient 

RepeatTor  next  medium 
Index  to  point  cross  sections  for 
each  multigroup  boundary 


3*NMIX 

NGP 

NGG 

NELEM*NC0EF 
if  IXTAPE  >  0 

NGP 

NGG 

NTG 

NTG 

NTG 

NGP 

NTG 

NTG 

NGG 

NGG*  NGP 
NGP* NGG 

NDSNGP 

NDSNGG 

NDSNGP*  NSCT 

NDSNGG*  NSCT 

NDSNGP*  NSCT 

NDSNGG*  NSCT 
ISP0RG 

NDSNGP 

NDSNGG 

INFP0G*(NC0EF-2) 

INF  P0G*(NC0EF- 1 ) 
(IGQPT+1)*NMED  if 
I06RT>O;  otherwise 
zero 


i\ 


I 
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"TS'S* 


’TOs^bwuwwmotuiw 


Location 

Information 

Size 

NXSECT(l) 

Point  cross 
sections  for 
Medium  1 

total 

scattering 

v'Li 

NPT(1)*(NEGPS-1)* 

NXPM 

NXSECT(J) 

* 

Point  cross 

total 

XT  T  ACT 

sections  for 
Medium  J 

scattering 
v  E  j 

NPT(J)*(NEGPS-1)* 

NXPM 

In  Jaiio  1 

Fig.  16.  Layout  of  permanent  cross  sections  in  blank  common. 
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■HSJTWW*' 


Location 


Variable 


TN0  =  NSIGL 


NN0  +  14 


IG 

U 

V 

w 

X 


Y 

z 

WATE 

AGE 

BLZNT 

NAME 

NAMEX 

NMED 

NREG 


Variables  are  from 
NUTR0N  common, 
see  Table  VI  for 
definitions. 


Repeat  for  Particle  2 


NN0  +  14*NM0ST  =  NLAST 

=  NFISBN  IF  MFISTP  >  0 


Fig.  17.  Layout  of  particle  bank  in  blank  common. 


.ajp  i  wmi  i WUI W.l I-  1  P,Uf*  *J*W*M*V'W^^*»*l*l 


Location 


Variable 


NFISBN 


NFISBN  +  7 


X 

Y 

Z 

WATEF 

AGE 

IG 

NAMEX 


Repeat  for  particle  2 


( 

S 


» 


NFISBN  +  7*NM0ST  =  NLAST 


Fig.  18.  Layout  of  fission  bank  in  blank  common. 


Location  Labels 

Mnemonic 
Variable  Name 

Length 

NLAST*  +  1 '  ' 

“  ■■ " 

L0CRSP 

LABELS 

20*NRESP 

RESP 

NRESP*NMTG 

L0CXD 

EXTR 

NEX*NMTG 

XD 

YD 

ZD 

RAD 

T0 

FACT 

6*ND 

L0CIB 

EXTR 

NEXND*ND 

IB 

EP 

3*NE 

L0CC0 

DELE 

L0CT 

C0S 

NA 

L0CUD 

T 

DELT 

2*ND*NT 

UD 

SUD 

3*ND*NRESP 

L0CSD 

SUD2 

SD 

SSD 

SSD2 

3*ND*NRESP 

L0COE 

SUD  &  SSD  Units 

20 

NOTE:  NLAST  must  not  be  changed  after  SC0RIN  is  called. 
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Location  Labels 


Mnemonic. 
Variable  Name 


Length 


QE 

SQE 

SQE 


3*NE*ND 


L0CQT. 


L0CQTE. 


L0CQAE 


LMAX. 


SQE  Units 


QT 

SQT 

SQT2 


SQT  Units 


QTE 

SQTE 

SQTE2 


SQTE  Units 


QAE 

SQAE 

SQAE2 


SQAE  Units 


20 


3  *NT*ND*NRESP 


20 


3*NT*NE*ND 


20 


3*NA*NE*ND 


20 


Fig.  19.  Layout  of  user  area  (analysis)  in  blank  common. 
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+  (NREG-1)*MAXGP  +  IG 


TABLE  ID.  (Cont'd. ) 
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TABLE  IV.  (Cont’d. ) 


TABLE  V. 


Definition  of  Variables  in  AP0LL0  Common 

Variable  Definition 

AGSTRT  Input  starting  age  of  source  particle. 

DDF  Starting  particle  weight  as  determined  in  S0RIN. 

DEADWT(5)  The  summed  weights  of  the  particles  at  death.  The 

four  deaths  are:  Russian  roulette,  escape,  energy, 
and  age  limit.  DEADWT(5)  is  unused. 

ETA  Mean -free -path  between  collisions. 

ET/.TH  Distance  in  cm  to  the  next  collision  if  the  particle 

does  not  encounter  a  change  in  total  cross  section. 

ETAUSD  Flight  path  in  m.  f.  p.  that  has  been  used  since  the 

last  event. 

UINP,  VINP,  WINP  Input  direction  cosines  for  source  particle. 

WTSTRT  Input  starting  weight. 

XSTRT,  YSTRT, 

ZSTRT  Input  starting  coordinates  for  source  particle. 

TCUT  Age  limit  at  which  particles  are  retired. 

XTRA(IO)  Not  used. 

10,  II  Output  and  input  logical  units. 

MEDIA  Number  of  media  for  which  there  are  cross  sections. 

IADJM  Switch  indicating  an  adjoint  problem  if  >  0. 

ISBIAS  Switch  indicating  that  source  energy  distribution  is 

to  be  biased  if  >  0. 

IS0UR  Input  source  energy  group  if  >  0;  otherwise,  S0RIN 

is  called  to  read  input  spectrum. 

ITERS  Number  of  batches  still  to  be  processed  in  the  run. 

I  TIME  Not.  used. 

ITSTR  Switch  indicating  that  secondary  fissions  are  to  be 

the  source  for  the  next  batch  if  >  0. 
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TABLE  V.  (Cont'd.) 


Variable 


Definition 


L0CWTS 

L0CFWL 

L0CEPR 

L0CNSC 

L0CFSN 

MAXGP 


MA^IM 


MEDALB 

MGPREG 

MXREG 

NALB 

NDEAD(5) 

NEWNM 

NGE0M 

NGPQT1* 

NGPQT2* 

NGPQT3* 


Starting  location  in  blank  common  of  the  weight  standards 
and  other  arrays  MGPREG  long . 

Starting  location  in  blank  common  of  the  fission  weights. 

Starting  location  in  blank  common  of  the  energy -biasing 
parameters. 

Starting  location  in  blank  common  of  the  scattering 
counters. 

Starting  location  in  blank  common  of  the  fission  and 
gamma-generation  probabilities  for  ea-.h  medium  and 
group. 

Maximum  number  of  energy  groups  for  which  there  are 
weight  standards  or  path-length  stretching  parameters. 

The  elapsed  clock  time  at  which  the  problem  is  termin¬ 
ated. 

Medium  number  for  the  albedo  medium. 

Product  of  number  of  weight  standard  groups  (MAXGP) 
and  regions  (MXREG). 

Maximum  number  of  regions  in  the  system. 

An  index  indicating  that  an  albedo  scattering  has  occurred 
if  >  0. 

Number  of  deaths  of  each  type  (see  DEADWT). 

Name  of  the  last  particle  in  the  bank. 

Location  of  first  cell  of  geometry  data  storage  in  blank 
commc 

The  lowest  energy  group  (largest  group  number)  for 
which  primary  particles  are  to  be  followed. 

The  number  of  primary  particle  groups. 

The  lowest  energy  group  (largest  group  number)  for 
which  any  particle  is  to  be  followed. 


* 

See  Appendix  A  for  diagram  of  energy  group  structure. 

77 


TABLE  V.  (Cont'd. ) 

Variable  Definition 

NGPQTG  Number  of  energy  groups  of  secondary  particles  to  be 

followed. 

* 

NGPQTN  Number  of  energy  groups  of  primary  particles  to  be 

followed. 

NITS  Number  of  batches  per  run. 

NKCALC  The  first  batch  to  be  used  for  a  k  calculation.  If  0,  k 

is  not  calculated. 

NK3LL  An  index  to  indicate  that  Russian  roulette  is  to  be 

played  if  >  0. 

NLAST  The  last  cell  in  blank  common  that  was  used  by  the 

cross-section  storage  or  is  set  aside  for  banking. 

NMEM  The  location  of  the  next  particle  in  the  bank  to  be  pro  - 

cessed. 

* 

NMGP  The  number  of  primary  particle  groups  for  which  there 

are  cross  sections. 

NM0ST  The  maximum  number  of  particles  that  the  bank  can  hold. 

* 

NMTG  The  total  number  of  energy  groups  (both  primary  and 

secondary)  for  which  there  are  cross  sections. 

N0LEAK  An  index  which  indicates  that  nonleakage  path-length 

selection  is  to  be  used  if  >0. 

N0RMF  An  index  to  indicate  that  the  fission  parameters  are  to 

be  renormalized  if  >0. 

NPAST  An  index  to  indicate  that  the  exponential  transform  is 

to  be  used  if  >  0. 

NPSCL(13)  An  array  of  counters  of  events  for  each  batch: 

1.  sources  generated 

2.  splittings  occurring 

3.  fissions  occurring 

4.  gamma  rays  generated 

_ 

See  Appendix  A  for  diagram  ol  energy  group  structure. 
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TABLE  V.  (Cont’d. ) 

Variable  Definition 


NQUIT 

NSIGL 

NS0UR 

NSPLT 

NXTRA(l-5) 

NXTRA(6) 

NXTRA(7-10) 


5.  real  collisions 

6.  albedo  scatterings 

7.  boundary  crossings 

8.  escapes 

9.  energy  cutoffs 

10.  time  cutoffs 

11.  Russian  roulette  kills 

12.  Russian  roulette  survivors 

13.  gamma  rays  not  generated  because  bank  was  full. 
Number  of  runs  still  to  be  processed. 

Starting  location  of  the  bank  in  blank  common. 

An  index  input  to  indicate  that  fissions  are  to  be  the 
source  for  future  batches. 

An  index  to  indicate  that  splitting  is  to  be  considered 
if  >0. 

Used  to  separate  INPUT1  and  INPUT2. 

NGPREG  calculated  in  M0RSE  and  used  in  G0MST  as 
an  index  for  track  length  estimation. 

Not  used. 


TABLE  VI. 


Definition  of  Variables  in  NUTR0N  Common 


Variable 

Definition 

NAME 

Particle's  first  name. 

NAMEX 

Particle's  family  name.  (Note  that  particles  do  not 
marry. ) 

IG 

Current  energy  group  index. 

IG0 

Previous  energy  group  index. 

NMED 

Medium  number  at  current  location. 

MED0LD 

Medium  number  at  previous  location. 

NREG 

Region  number  at  current  location. 

U,  V,  W 

Current  direction  cosines. 

U0LD,  V0LD, 
W0LD 

Previous  direction  cosines. 

X,  Y,  Z 

Current  location. 

WATE 

Current  weight. 

0LDWT 

Weight  at  previous  collision. 

WTBC 

Weight  just  before  current  collision. 

BLZNT 

Current  block  and  zone  number  (packed)  for  05R  geo¬ 
metry,  and  current  zone  number  IR  for  combinatorial 
geometry. 

BLZ0N 

Previuus  block  and  zone  number  (packed)  for  05R  geo¬ 
metry,  and  previous  zone  number  IR  for  combinatorial 
geometry. 

AGE 

Current  age. 

0LDAGE 

Previous  age. 

TABLE  Vn. 


Definitions  of  Variables  in  USER  Common 
Variable  Definition 


AGSTRT 

WTSTRT 

XSTRT 

YSTRT 

ZSTRT 

DFF 

EB0TN 

EB0TG 

TCUT 

10 

II 

IADJM 

NGPQT1 

NGPQT 

NGPQT3 

NGPQTG 

NGPQTN 

NITS 

NLAST 

NLEFT 


Initial  chronological  age  to  be  assigned  to  source 
particles. 

Initial  weight  to  be  assigned  to  source  particles. 

Initial  x  position  to  be  assigned  to  source  particles. 

Initial  y  position  to  be  assigned  to  source  particles. 

Initial  z  position  to  be  assigned  to  source  particles. 

Normalization  for  adjoint  problems. 

Lower  energy  boundary  (eV)  of  last  neutron  group 
(group  NMGP). 

Lower  energy  boundary  (eV)  of  last  gamma -ray  group 
(group  NMTG). 

Chronological  age  limit. 

Logical  unit  for  output. 

Logical  unit  fcr  input. 

Adjoint  switch  (>  0  for  adjoint  problem). 

Problem -dependent  energy  group  limits  (see  Appendix 
A). 

Group  number  of  lowest  energy  gamma-ray  group  to  be 
treated. 

Group  number  of  lowest  energy  neutron  group  to  be 
treated. 

Number  of  batches  to  be  run. 

Last  cell  in  blank  common  used  by  random  walk  package. 
Number  of  cells  in  blank  common  available  to  user. 


NMGP  Number  of  primary  particle  energy  groups. 

NMTG  Total  number  of  energy  groups. 
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TABLE  VIII. 

Definitions  of  Variables  in  Common  G0ML0C 


Variable 


Definition 


KMA  Starting  location  for  the  array  MA  containing  integer 

data  for  each  code  zone. 

KFPD  Starting  location  for  the  array  FPD  containing  real  data 

for  each  code  zone. 

KLCR  Starting  location  for  the  array  L0CREG(I)  which  contains 

the  starting  location  in  the  MA  array  for  the  Ith  code 
zone  data. 

KNBD  Starting  location  for  the  array  NUMB0D(I)  which  con¬ 

tains  the  number  of  bodies  for  the  I*h  code  zone. 

KI0R  Starting  location  for  the  array  IR0R(I)  which  contains 

the  index  of  the  corresponding  input  zone  for  the  Ith 
code  zone. 

KKIZ  Starting  location  for  the  array  MIRZ(I)  which  contains 

the  index  of  the  M0RSE  region  corresponding  to  the  Ith 
geometry  input  zone. 

KRCZ  Starting  location  for  the  array  MRCZ(I)  which  contains 

the  index  of  the  M0RSE  region  corresponding  to  the 
geometry  code  zone. 

KMIZ  Starting  location  for  the  array  MMIZ(I)  which  contains 

the  index  of  the  M0RSE  media  corresponding  to  the  I^1 
geometry  input  zone. 

KMCZ  Starting  location  for  the  array  MMCZ(I)  which  contains 

the  index  of  the  M0RSE  media  corresponding  to  the  I**1 
code  zone. 

KKR1  Starting  location  for  the  array  KR1(L)  contains  the  first 

code  zone  which  was  made  from  the  input  zone. 

KKR2  Starting  location  for  the  array  KR2(L)  contains  the  last 

code  zone  which  was  made  from  the  input  zone. 

KNSR  Starting  location  for  array  NS0R  which  contains  the  code 

zones  in  which  source  particles  have  been  found. 
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TABLE  Vm.  (Cont’d.) 


Variable 

Definition 

KV0L 

Starting  location  for  the  array  VN0R(I)  which  contains 
the  volume  for  M0RSE  region  (I). 

NADD 

Starting  location  for  the  geometry  data  length  and 
changed  in  J0MIN  to  the  total  number  of  words  re¬ 
quired  for  geometry  data. 

LDATA 

Length  of  the  integer  data  in  the  MA  array  excluding 
the  words  set  aside  for  zone  search  information. 

LTMA 

Total  length  of  the  MA  array. 

LFPD 

Length  of  the  FPD  array 

NUMR 

Number  of  code  produced  zones. 

IRTRJ 

Number  of  input  zones. 

NUMB 

Number  of  bodies. 

NIR 

Number  of  M0RSE  geometry  regions. 

y*  1 M  Wf’ly^Jr 
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TABLE  IX. 


Definitions  of  Variables  in  Common  L0CSIG 


Variable 

Definition 

ISTART 

Starting  location  for  the  total  cross-section  vector  for 
the  first  medium. 

ISCC0G 

Starting  location  i„  .  the  scattering  cross-section  veccor 
for  the  first  medium. 

INAB0G 

Starting  location  for  the  non-absorption  vector  for  the 
first  medium. 

IGAB0G 

Starting  location  for  the  gamma-ray  production  vector 
for  the  first  medium. 

IFP0RG 

Starting  location  for  v  E  f  ,  the  fission  neutron  production 
vector  for  the  first  medium. 

IFNGP 

Starting  location  for  the  primary-secondary  transfer 
probability  matrix. 

IFSP0G 

Starting  location  of  the  primary  downscatter  probability 
matrix. 

IDSG0G 

Starting  location  of  the  secondary  downscatter  prob¬ 
ability  matrix. 

IPRBNG 

Starting  location  of  the  primary  scattering  angle  prob¬ 
ability  matrix. 

IPRBGG 

Starting  location  of  the  secondary  scattering  angle  prob¬ 
ability  matrix. 

ISCANG 

Starting  location  of  the  primary  scattering  angle  matrix. 

ISCAGG 

Starting  location  of  the  secondary  scattering  angle 
matrix. 

ISP0RG 

Size  of  storage  needed  for  each  medium,  not  including 
Legendre  coefficients. 

vymrmxtv* *?<r  r^rip^i  ■•»  wj»  • w ■u^y^ ■.! m  f  -|!M« W+, VS v*?J  w.'v'  j>  >n  »j. 


1 F  *»  -WV. 


TABLE  IX.  (Cont'd. ) 

Variable  Definition 


ISP0RT* 

INPBUF 

ISIG0G 

INFP0G* 

IABS0G 

IT0TSG* 

NGP 

NDS 

NGG 

NDSG 

INGP 

INDS 

NMED 

NELEM 


Starting  location  for  temporary  storage  of  downscatter 
matrix. 

Starting  location  for  temporary  storage  of  the  P^  table. 

Starting  location  for  temporary  storage  of  total  cross 
section  for  element  1. 

Starting  location  for  temporary  storage  of  ul»f  for 
element  1. 

Starting  location  for  temporary  storage  of  downscatter 
matrix  for  P^  coefficients  (primary  groups,  element  1). 

Totai  storage  required  for  temporary  storage. 

The  number  of  primary  groups  to  be  treated. 

Number  of  downscatters  for  NGP  (usually  equal  to  NCP). 

Number  of  secondary  groups  to  be  treated. 

Number  of  downscatters  for  NGG  (usually  equal  to  NGG). 

Number  of  groups  for  which  cross  sections  are  to  be 
inp  \t. 

Number  of  downscatters  for  the  INGP  groups. 

Number  of  media  for  which  cross  sections  are  to  be 
stored --should  be  same  as  MEDIA  as  read  on  Card  B 
of  M0RSE  input. 

Number  of  elements  for  which  cross  sections  are  to  be 
read. 


If  Legendre  coefficients  are  to  be  restored,  then: 

INFP0G  -  Redefined  by  JNPUT  as  number  of  locations  required 
for  each  coefficient  (both  primary  and  secondary). 

IT0TSG  -  Redefined  by  JNPUT  as  total  storage  required  for  all 
coefficients  for  each  medium. 

ISP0RT  -  Redefined  by  JNPUT  as  starting  location  of  Pj  co¬ 
efficient  for  primary  groups  for  medium  1. 


85 


TABLE  IX.  (Cont'd.) 


Variable 


Definition 


NMIX 

NC0EF 

NSCT 

NTS 

NTG 

NDSNGP 

NDSNGG 

IADJ 

NME 

L0C 

INGS 

INSG 

II,  10 
KKK 

IXTAPE 

IDEL 


Number  of  elements  times  density  operations  to  be  per 
formed. 


Number  of  coefficients,  including  Pq  . 

Number  of  discrete  angles  (usually  NC0EF/2Jntegraj). 

Number  of  downscatters  for  combined  primary  and 
secondary  groups  (usually  equal  to  NTG). 


Total  number  of  groups  (primary  +  secondary)  =  NGP 
f  NGG. 


The  number  of  locations  needed  for  tne  downscatter 
matrix  for  the  primary  particles. 

The  number  of  locstions  needed  for  the  downscatter 
matrix  for  the  seconc  \ry  particles. 

Same  as  IADJM. 


Indicator  for  stripping  gamma  rays  from  a  coupled 
neutron  gamma-ray  cross-sectior  set-~set  equal  to 
number  of  neutrons  groups  +  1. 

Same  us  L0CEPR. 

Starting  location  of  the  indices  fi  r  starting  location  of 
the  downscatter  vector  for  each  group  for  primary 
particles. 

Starting  locat;-  of  the  indices  for  starting  location  of 
tne  downscatter  vector  for  each  group  for  secondary 
particles. 

Input  and  output  logical  unit  numbers. 

A  running  index  of  the  number  of  cross  sections  that 
have  already  been  read  in  (used  in  checking  the  element 
numbers  obtained  from  tape). 


Logical  tape  number  of  the  multigroup  cross-section 
tape  if  >  0,  or  logical  tape  number  of  the  processed 
cross-section  tape  if  <  0. 

Starting  location  for  element  identifiers  which  determine 
the  element  cross  sections  to  be  read  from  tape. 
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TABLE  IX.  (Cont'd.) 


Variable  Definition 

ITEML  Amount  of  storage  for  primary  and  secondary  group 

downscatters  per  element. 

ITEMG  Starting  location  for  temporary  storage  of  downscatter 

matrix  for  Pl  coefficients  (secondary  groups)  for 
element  1. 

IRSG  Starting  location  of  the  mixing  parameters. 

IRDSG  Switch  to  print  the  cross  sections  and  to  test  the  card 

sequence  as  they  are  read  if  >  0  (test  card  sequence 
only  if  =  0,  and  does  neither  if  <0). 

ISTR  Switch  to  print  cross  sections  as  they  are  stored  if 

>0. 

IPRIN  Switch  to  print  angles  and  probabilities  if  >0. 

IFMU  Switch  to  print  intermediate  results  of  n's  calculation 

if  >  0. 

IM0M  Switch  to  print  moments  of  angular  distribution  if  >  0. 

IDTF  Switch  to  signal  that  input  format  is  DTF-IV  format  if 

>0;  otherwise,  ANISN  format  is  assumed. 

ISTAT  Flag  to  restore  Legendre  coefficients  for  next -flight 

estimates  if  >0. 

IPUN  Switch  to  print  results  of  bad  Legendre  coefficients  if 

>0. 

NUS  Number  of  groups  of  upscatter. 

NGN  Not  use  A  presently. 

IHT  Location  of  S-p  in  the  cross-section  table,  currently 

set  to  3  in  XSEC. 

INUS  Starting  location  for  the  upscattering  vector  for  the 

first  medium. 

INUSN  Number  of  impossible  first  moments  found  in  the  cross 

sections. 

INGN  Starting  location  for  the  photoneutron  production  vector 

for  the  first  medium  (not  used  at  present--same  as 
INUS). 
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TABLE  IX.  (Cont'd.) 

Definition 

Starting  location  for  the  secondary -primary  transfer 
probability  matrix  (not  used  at  present --same  as  INUS). 

Starting  location  for  the  array  of  the  number  of  down- 
scatter  groups  for  each  primary  group. 

Starting  location  for  the  array  of  the  number  of  down- 
scatter  groups  for  each  secondary  group. 


Variable 

INGNP 

INNN 

IGGG 


^T-TT"1  T'*  ?TTr 


m  ■  i  Wpp^i  wgt*  w*i ,n 


rryy 1  mu J.W  nflmuM'Hm }  w»w W, 

/ 


TABLE  X. 

Definitions  of  Variables  in  Common  GTSC1 


Variable 


Definition 


NPT(16) 


NXSECT(17) 

NDSGP 
II,  10 
I06RT 
IGQPT 

NEGPS 

ESPD(IOO) 

NNIC 

NXPM 


NINC 


Number  of  points  per  supergroup  for  each  medium. 
(Only  16  media  are  allowed. )  Redefined  at  end  of 
GTSCT  as  the  total  storage  required  for  each  cross 
section  for  each  medium. 

Starting  location  of  point  cross  sections  for  each 
medium. 

Total  storage  required  per  supergroup  for  all  media. 

Input  and  output  logical  unit  numbers. 

Logical  tape  number  of  06R  point  cross-section  tape. 

Last  M0RSE  multigroup  for  which  point  cross  sections 
will  be  used  NMGF'. 

Number  of  supergroup  boundaries. 

Energy  boundaries  (in  eV)  of  the  supergroups  (only  100 
supergroups  are  allowed). 

Starting  location  of  the  point  cross  sections. 

Number  of  cross  sections  per  medium 

=  1  if  total  only 

=  2  if  total  +  scattering 

=  3  if  total,  scattering,  and  y*fission. 

Index  for  locating  point  cross  sections.  Used  as  flag  in 
FISGEN  to  determine  if  transport  process  has  begun. 
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TABLE  XI. 

Definitions  of  Variables  in  Common  PDET 


Variable 

Definition 

ND 

Number  of  detectors  (2  1). 

NNE 

Number  of  neutron  energy  bins. 

NE 

Total  number  of  energy  bins. 

NT 

Number  of  time  bins. 

NA 

Number  of  angle  bins. 

NRESP 

Number  of  energy  response  functions  1). 

NEX 

Number  of  extra  arrays  of  length  NMTG  to  be  set  aside 

NEXND 

Number  of  extra  arrays  of  length  ND  to  be  set  aside. 

NEND 

NE*ND. 

NDNR 

ND*NRESP. 

NTNR 

NT*NRESP. 

NTNE 

NT*NE. 

NANE 

NA*NE. 

NTNDNR 

NT*ND*NRESP. 

NTNEND 

NT*NE*ND. 

NANEND 

NA*NE*ND. 

L0CRSP 

Location  of  cell  zero  of  response  functions. 

L0CXD 

Location  of  cell  zero  of  detector  positions. 

L0CIB 

Location  of  cell  zero  of  energy  bins. 

L0CC0 

Location  of  cell  zero  of  angle  (cosine)  bins. 

L0CT 

Location  of  cell  zero  of  time  bins. 

L0CUD 

Location  of  cell  zero  of  array  UD. 

L0CSD 

Location  of  cell  zero  of  array  SD. 

L0CQE 

Location  of  cell  zero  of  array  QE. 

L0CQT 

Location  of  cell  zero  of  array  QT. 
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TABLE  XI.  (Cont'd.) 

Variable  Definition 

L0CQTE  Location  of  cell  zero  of  array  QTE. 

L0CQAE  Location  of  cell  zero  of  array  QAE. 

LMAX  Last  cell  used  in  blank  common. 

EFIRST  Upper  energy  limit  of  first  energy  bin. 

EGT0P  Upper  energy  limit  of  first  gamma-ray  bin. 


TABLE  Xn. 

Definitions  of  Variables  in  Common  FISBNK 

Variable  Definition 

MFISTP  Index  to  indicate  that  fissions  are  to  be  treated  if  >  0. 

NFISBM  Location  in  blank  common  of  cell  zero  of  the  fission 

bank. 

NFISH  Number  of  fissions  produced  in  the  previous  batch 

(generation). 

FT0TL  Total  fission  weight  from  all  collisions  in  this  batch 

(generation). 

FWATE  Total  weight  of  fission  neutrons  stored  in  bank  in  this 

batch  (generation). 

Weight  of  fission  neutron  stored  in  bank. 


WATEF 


TABLE  Xm. 


definitions  of  Variables  in  Common  GE0MC 
as  found  in  05R  Geometry 

Variable  Definition 

X2,  Y2,  Z2  Coordinates  at  tentative  end -of -flight  or  if  the  trajectory 
is  in  an  internal  void;  X2,  Y2,  Z2  are  the  direction 
cosines  of  the  trajectory. 

XI,  Yl,  Zl  Starting  coordinates  for  the  particle. 

ETA  Number  of  mean  free  paths  to  be  traversed  if  flight  goes 

to  X2,  Y2,  Z2. 

ETA  USD  Number  of  mean  free  paths  actually  traversed  after  the 

call  to  GE0M. 

IBLZ  An  index  to  the  medium  number  for  the  special  geometry 

packages.  For  GENERAL  GE0M,  IBLZ  is  a  packed 
work  giving  the  block  and  zone  of  the  end  of  flight. 

IBZN  A  dummy  variable. 

MARK  A  flag  set  by  GE0M  indicating  the  results  of  the  trajec¬ 

tory  calculation 

=  1  for  completed  flight 
=  0  for  boundary  crossing 
=  -1  for  escape 

=  -2  for  entering  an  internal  void. 

NMED  Medium  number  at  end  of  flight  or  of  medium  about  to 

be  entered  at  a  boundary  crossing. 

Region  number  at  end  of  flight;  not  set  at  boundary 
crossings. 


NREG 
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TABLE  XIV. 

Definitions  of  Variables  in  Common  PAREM 
as  found  in  Combinatorial  Geometry 

Variable  Definition 


XB(3) 

WB(3) 

IR 

WP(3) 

XP(31 

IDBG 

IRPRIM 

NASC 

LSURF 

NB0 

LRI 

LR0 

RIN 

R0UT 

KL00P 

L00P 

ITYPE 


Coordinates  of  the  starting  point  of  the  present  path. 

Direction  cosines  of  particle  trajectory.  Equal  to  U, 

V,  and  W. 

Combinatorial  zone  of  present  particle  position. 

Temporary  storage  of  WB(3). 

Temporary  storage  of  XB(3). 

Set  non-zero  to  initialize  a  debug  printout. 

Next  region  to  be  entered  after  a  call  of  Gl. 

Body  number  of  last  calculated  intersection.  Set  negative 
to  indicate  source  or  collision  point  not  on  a  body  sur¬ 
face. 

Surface  of  body  NASC  where  intersection  occurred. 
Positive  if  particle  is  entering  the  body  and  negative 
when  exiting. 

Body  number  and  a  sign  used  to  define  zones.  Input  in 
zone  description  as  positive  when  zone  is  contained  in 
body  and  as  negative  if  zone  is  outside  body. 

Entry  surface  calculated  in  GG. 

Exit  surface  calculated  in  GG. 

Distance  to  entry  calculated  in  GG. 

Distance  to  exit  calculated  in  GG. 

Trajectory  index  of  present  path  incremented  in  Gl. 

Index  of  last  trajectory  calculated  for  body  NB0 .  If 
L00P  is  equal  to  KL00P,  GG  returns  immediately  with 
old  values  of  RIN,  R0UT,  LRI,  and  LR0. 

Body  type  of  body  NB0  (indicates  B0X,  SPH,  etc. ). 
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TABLE  XIV.  (Cont’d. ) 


r 

• 

Variable 

Definition 

L 

PINF 

Machine  infinity. 

N0A 

Not  used. 

f 

DIST 

Distance  from  XB(3)  to  present  point. 
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TABLE  XV. 

* 

Definitions  of  Variables  in  Common  0RGI 

Variable  Definition 

M.L.  jC  Set  1  in  G1  if  trajectory  end  point  is  reached  before 

next  intersection.  Otherwise  set  to  0. 

DISTO  Distance  from  point  XB(3)  to  next  scattering  point. 

Used  in  G1  to  avoid  calculating  the  next  zone  if  a 
scattering  event  occurs  before  the  intersection. 

NMEDG  Zone  number  IR  from  a  L00KZ  call.  Stored  in  BLZNT 

by  MS0UR. 

_ 

Note:  Variabl  names  are  not  the  same  in  all  routines. 


96 


TABLE  XVI. 


Definitions  of  Random  Walk  Variables  in  Blank  Common 


Mnemonic 

Variable  Definition 

Name 


ENER(IG) 

VEL(IG) 


FS(IG) 


BFS(IG) 

WTHI(IG, 

NREG) 

WTL0(IG, 

NREG) 

WTAV(IG, 

NREG) 

PATH(IG, 

NREG) 

NSPL(IG, 

NREG) 

WSPL(IG, 

NREG) 

N0SP(IG, 

NREG) 

WN0S{IG, 

NREG) 

RRKL(IG, 

NREG) 


Upper  energy  boundary  of  group  IG  (in  eV). 

Velocity  corresponding  to  the  mean  energy  for  neutron 
groups  and  the  speed  of  light  for  gamma-ray  groups 
(in  cm/sec). 

Unbiased  source  spectrum — unnormalized  fraction  of 
source  particles  in  each  energy  group --transformed  to 
c.d.  f.  by  S0RIN. 

Biased  source  spectrum — relative  importance  of  each 
energy  group— transformed  to  biased  c.  d.  f.  by  S0RIN. 

Weight  above  which  splitting  is  performed  (vs.  group 
and  region). 

Weight  below  which  Russian  roulette  is  performed  (vs. 
group  and  region). 

Weight  tc  be  assigned  Russian  roulette  survivors  (vs. 
group  and  region). 

Exponential  transform  parameters  (vs.  group  and  region). 
Splitting  counter  (vs.  group  and  region). 

Weight  equivalent  to  NSPL. 

Counter  for  full  bank  when  splitting  was  requested  (vs. 
group  and  region). 

Weight  equivalent  to  N0SP. 

Russian  roulette  death  counter  (vs.  group  and  region). 
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TABLE  XVI.  (Cont'd.) 


Mnemonic 

Variable  Definition 

Name 


WRKL(IG, 

NREG)  Weight  equivalent  to  RRKL. 

RRSU(IG,  Russian  roulette  survival  counter  (vs.  group  and 

NREG)  region). 

WRSU(IG, 

NREG)  Weight  equivalent  to  RRSU. 

INIWHI  (IG, 

NREG)  Initial  values  of  WTHI  array. 

INIWL0(IG, 

NREG)  Initial  values  of  WTL0  array. 

INIWAV(IG, 

NREG)  Initial  values  of  WTAV  array. 

FWL0(NREG)  Weights  to  be  assigned  to  fission  daughters  (vs. 

region). 

INIFL0(NREG)  Initial  values  of  FWL0. 

GWL0(IG,  Weights  to  be  assigned  to  secondary  particles  (vs. 

NREG)  group  and  region). 

EPRB(IG,  Relative  importance  of  energy  groups  after  scatter- 

NREG)  ing  (vs.  group  and  region). 

NSCT(IG, 

NREG)  Number  of  real  scatterings  (vs.  group  and  region). 

WSCT(IG, 

NREG)  Weight  equivalent  to  NSCT. 

NALB(IG,  Number  of  albedo  scatterings  (vs.  group  and 

NREG)  region). 

WALB(IG, 

NREG)  Weight  equivalent  to  NALB. 

NFIZ(IG, 

NREG)  Number  of  fissions  (vs.  group  and  region). 
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lABLEXVI.  (Cont'd. ) 

Mnemonic  * 

Variable  Definition 

Name 


WFIZ(IG, 

NREG)  Weight  equivalent  to  NFIZ. 

NGAM(IG,  ,  ' 

NREG)  Number  of  secondary  productions  (vs.  group  and  region). 

WGAM(IG, 

,NREG)  ,  Weight  equivalent  to  NGAM.  i 

NSCA(IMED)  Scattering  counter  (vs.  cross-section  medium). 

FISH(IG,  Probability  of  generating  fission  neutron  (vs.  group 

IMED)  and  medium). 

FSE(IG,  Source  spectrum  for  fis$ion -induced  neutrons  for  each 

IM]ED)  group— input  as  frequency  of  group  IG. 

GMGN(IG,  Probability  of  generating  secondary  particle  (vs.  group 

IMED)  and  medium). 


THE  PICTURE  PROGRAM 


IV. 

\ 

This  program  is  used  to  obtain  printer  plots  of  two-dimensional 
slices  through  the  three-dimensional  geometry  described  by  the  combi¬ 
natorial  geometry  (CG)  package.  The  05R  geometry  version^°|  has 
been  modified  for  use  with  CG  and  some  minor  changes  to  the  program 
have  been  made.  It  is  now  possible  to  obtain  a  picture  the  full  width  of 
a  computer  page  and  the  logic  for  determining  the  medium  and/or  region 
number  has  been  changed  so  that  a  test  of  the  geometry  ray  tracing  cap¬ 
ability  is  also  made.  That  is.  instead  of  just  testing  the  routine  L00KZ, 
this  modification  forces  the  geometry  package  to  execute  the  routines 
and  logic  used  during  the  random  walk  process.  This  program  has  been 
very  useful  in  debugging  geometry  input.  It  is  suggested  tjiat  any  com¬ 
plex  geometry  be  tried  in  this  program  before  using  it  in  the  M0RSE 
code.  A  sample  problem  is  provided.  Besides  the  geometry  package, 
a  main  program,  PICTURE,  and  two  routines,  MESH  and  PRINT,  are  re¬ 
quired.  A  discussion  of  these  routines  is  included  and  a  definition  of 
the  variables  in  common ,  PICT,  is  given  in  Table  XVII. 


Preceding  page  blank 
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4. 1  Program  PICTURE 

PICTURE  is  the  executive  routine  for  the  PICTURE  program 
and  reads  in  the  input,  calculates  the  coordinates  of  the  picture  to  be 
plotted  and  controls  the  calls  to  other  routines.  There  are  several  dif¬ 
ferent  ways  in  which  a  two-dimensional  slice  through  the  geometry  may 
be  obtained.  These  different  options  are  discussed  in  Appendix  B  with 
discussions  of  the  input.  The  characters  to  be  printed  for  corresponding 
media  or  regions  may  be  changed  by  altering  the  values  in  ATABLE  as 
given  in  the  data  statement. 

Subroutines  called:  J0MIN,  PRINT 
Commons  required:  PICT 

Variables  required:  Several  input  cards  are  read. 

Variables  changed:  All  variables  in  common  PICT. 

Significant  internal  variables: 


INT 

input  logical  unit, 

I0T 

output  logical  unit, 

NADD(l) 

first  location  in  blank  common  for  storage 
of  geometry  data,  . 

NCK 

flag  to  indicate  which  set  of  input  options 
was  used  to  define  the  two-dimensional 
slice. 

Program  PICTURE 


(£> 

®- 


Set  input/output  devices  to  5  ancTcT] 


I/£):  Symbol  option. 


Read  symbol  table  if  required 


Call  J0MIN  to  input  geometry  data. 


1/0:  Options,  number  of  pictures  for  this 
geometry,  media  or  region  title. 


T/0:  picture  coordinates  upper  left  and 


lower  right  corner. 

- 


1/0:  direction  numbers  for  U  and  V  axes^ 


ft/07  interval  number  and 


size. 


T 

£ 


- —  - < 

No 

,KCK=1 

Use_input  DELU  or  DELV  values  and  NU  Nvl 

|C.  !!  PRINT  ] 

.  _ 1 

-  \  is  there  another  geometry  input’} 

Calculate 
DELU  given 
NU  or  NV 
given  DELU 
repeat  for  V 


4.2 


Subroutine  PRINT  (KXX,  KYY,  ATABLE) 

This  routine  controls  the  printing  of  the  picture.  First  it  is 
determined  if  the  picture  is  to  be  more  than  one  "page”  (130  columns) 
wide  and  then  for  each  line  of  the  picture  on  the  first  page  the  infor¬ 
mation  to  be  printed  is  determined  by  calling  MESH  and  printed  with  a 
130A1  format.  The  next  lines  are  then  calculated  and  printed  until  that 
"page”  is  finished.  Note  that  a  page  refers  to  width  not  length.  Thus, 
as  much  detail  may  be  obtained  as  necessary  in  both  directions  by  piec¬ 
ing  together  the  output. 

Called  from:  PICTURE 

Subroutine  called:  MESH 
Commons  required:  PICT 
Variables  required: 

KXX  -  number  of  intervals  in  the  U  direction 

(direction  of  paper  movement  through 
the  printer), 

KYY  -  number  of  intervals  in  the  V  direction 

(line), 

ATABLE  -  table  of  characters  to  be  printed. 

Significant  interval  variable: 

NPAGES  -  number  of  subpictures  required  to  cover 
the  width  of  the  total  picture, 

I0T  -  output  logical  unit, 

NV  -  number  of  characters  per  line  (characters/ 

page  width). 
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re*-.-  «,  If  *  ***  "A  ■-  '.’M  V,  M.Ai'.lU.  4-,; 


JU'IW 


*1 


Subroutine  PRINT 


RETURN 


4.  3  Subroutine  MESH  (XS,  YS,  ZS,  NV) 

Subroutine  MESH  is  used  by  the  PICTURE  package  to  set  up 
one  line  of  print  in  the  array  NST0R.  Both  for  efficiency  and  to  debug 
the  combinatorial  geometry  package,  this  version  has  been  modified  to 
work  exactly  like  particle  tracking.  L00KZ  is  first  called  to  deter¬ 
mine  the  zone  of  the  first  grid  point.  A  trajectory  to  the  last  grid  point 
is  then  initialized,  and  successive  calls  to  G1  "track”  a  particle  to  the 
last  point,  setting  the  region  of  each  grid  point  in  NST0R.  By  setting 
IRG  negative,  °ro,  or  positive,  either  NREG,  IR,  or  NMED  will  be 
stored  in  the  print  array  NST0R. 

Called  from:  PRINT 

Subroutines  called:  L00KZ,  G1 

Commons  required:  PICT,  G0ML0C,  PAREM,  0RGI,  blank. 

Variables  required: 

XS,  YS,  ZS  -  coordinates  of  first  grid  point, 

NV  -  number  of  grid  points, 

DELV  -  distance  between  grid  points, 

IRG  -  flag  to  print  NREG,  IR,  or  NMED  if  IRG 

is  negative,  zero,  or  positive.  . 

Variables  changed:  NST0R  -  print  array. 

Significant  internal  variables: 

J  -  grid  point  index, 

IST0R  -  value  to  be  stored  in  NST0R  between 

successive  boundary  crossings. 
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Subroutine  MESH 


[Call  L0PKZ  to  determine  IR  of  first  region. 


Store  NREG,  IR,  or  NMED  in  NSTQR  (1)  depending  on  IRG 


Initialize  combinatorial  geometry  trajectory  from] 
_ the  first  to  the  last  grid  p oi nt. _ 


Set  next  grid 


point  to  2, 


Call  Gl. 


Set  ISTOR  to  NREG,  IR,  or  NMED  depending  on  IRG 


From  DIST,  determine  last  grid  point  before 
boundary  crossing  or  trajectory  completion. 


Yes- 


Ir  last  grid  point  less  than  the  next  grid  point?  | 
(Distance  is  less  than  grid  resolution. ) _ / 


No 


1 

Loop  from  next  to 

last  .grid  points. 

Store  IST0R  in  NST0R  array. 

i _ 

_ 1 

Set  next  point  t( 
— - 

3  last  point  +1. 

flR  r  IR 

PRIM  1 

(Is  trajectory  complete?  ] - 

retur: 


TABLE  XVH. 


* 

Definition  of  Variables  in  Common  PICT 

Variable  Definition 

DELU  The  increment  in  geometry  units  between  lines  in  the 

picture  in  the  U  direction. 

DELV  The  increment  in  geometry  units  between  lines  in  the 

picture  in  the  V  direction. 

XO,  YO,  ZO  The  coordinates  in  geometry  units  defining  the  first 
point  (upper  left  hand  corner)  of  the  picture. 

XU,  YU,  ZU  The  length  (in  direction  U)  of  the  picture  in  geometry 

units. 

XV,  YV,  ZV  The  width  (in  direction  V)  of  the  picture  in  geometry 

units. 

NST0R(13O)  An  array  used  to  store  the  medium,  region,  or  zone 
number  for  one  line  of  the  picture. 

IRG  A  flag  indicating  that  region,  zone,  or  medium  geo¬ 

metry  parameter  should  be  printed  if  IRG  is  negative, 
zero,  or  positive,  respectively. 


V.  M0RSE  SAMPLE  PROBLEMS 

5. 1  Sample  Problem  Number  1 

(7) 

This  sample  problem  is  similar  to  that  used  previously 
which  calculates  neutrons  and  secondary  gamma  rays  in  infinite  homo¬ 
geneous  air.  The  source  is  a  12.  2-  to  15-MeV  neutron  source.  Two 
types  of  estimators  are  used  in  the  analysis  of  the  particle  histories. 
First,  a  boundary-crossing  estimator  is  used  for  five  spherical  shells. 
(This  part  of  the  problem  is  similar  to  the  M0RSE  sample  problem. ) 
Second,  a  collision  density  estimator  is  used  in  RELC0L  for  nine  de¬ 
tector  regions  (detectors  6-14).  By  choosing  the  spherical  shells  as 
the  midpoints  of  the  regions,  an  internal  check  is  provided.  Also,  by 
choosing  the  geometry  regions  to  correspond  to  the  collision  density 
detectors,  the  results  from  SAMB0  can  be  compared  with  results  from 
ENDRUN.  With  the  inclusion  of  the  combinatorial  geometry  package, 
an  additional  estimator  based  on  track  length/unit  volume  has  been 
added  with  results  also  output  in  ENDRUN.  Both  the  collision  density 
and  track  length  per  unit  volume  estimators  are  obtained  for  the  geo¬ 
metry  regions  without  s.  ecific  coding. 

Since  the  collision  density  estimator  in  RELC0L  scores 

2 

4trR.  <p  there  will  be  some  differences  between  these  estimates  and 

those  from  the  collision  ,.ensity  estimator  output  in  ENDRUN  which 
2 

scores  4ir  <R>  <p.  (The  symbols  <  >  denote  an  average  over  a  spatial 

region. ) 

Table  XVIII  illustrates  the  comparison  for  one  of  several 
sets  of  answers  as  calculated  by  the  various  estimators.  Descriptions 
of  the  user  routines  INSC0R,  RELC0L,  and  BDRYX  required  for  the 
sample  problem  are  given. 
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TABLE  XVm 


Comparison  of  Results  of  Sample  Problem  1 


Detector  No. 

Type  and  Location 
of  Detector 

Result 

Response  #1  Group  1  * 

4ffR2  fluence  fluence/eV 
top  group 

3 

Boundary  crossing 
at  450  meters. 

8. 19(-2) 

2.  92(-8)*A 
=  8. 19(-2) 

9 

Collision  density  for 
volume  between  400- 
500  meters. 

7.  59(-2) 

2.  71(-8)*A 
=  7.  59(-2) 

Region  4 

Volume  between 
400-500  meters 
(ENDRUN  results). 

Real 

Collision 
Region  4 

Top  Group 

1. 11(-18)*A*B 
=  7.  76(-2) 

* 

Volume  between 
400-500  meters 
(ENDRUN  results). 

Track 

Length 

Region  4 

Top  Group 

1. 18(-18)*A*B 
=  8.  36(-2) 

A  =  width  of  the  top  group  2.  8 (+6)  eV. 
B  =  4ff  {R:  ^  for  4^  region  =  2.  5(+10). 
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5.1.1 


Subroutine  INSC0R  -  Sample  Problem  Number  1 


This  user  routine  is  used  to  calculate  the  volume  of  each  of 
the  detector  regions  (detectors  6-14)  in  the  sample  problem.  The  re¬ 
ciprocal  of  the  volume  is  stored  in  the  array  in  blank  common  set  aside 
for  FACT,  a  detector -dependent  normalization.  Because  two  types  of 
estimators  are  used  in  the  problem,  normalization  values  of  1  are 
stored  for  the  boundary  detectors.  The  detector  number  NDC  for  the 
last  boundary  crossing  detector  (in  this  case,  5)  must  be  input. 

Called  from:  SC0RIN 

Subroutines  called:  None 

Commons  required:  PDET,  USER,  DETCUT,  blank. 

Variables  required: 

NDC  -  the  number  of  the  last  boundary  crossing 

detector.  Detector  radii  from  blank 
common. 

Variables  changed: 

Normalization  array  in  blank  common. 


Ill 


Subroutine  INSC0R 


?ART 


I:  NDC 


Yes 

I<ND 


'Are  there' 
more 

^detectors ?, 


{RETURN) 


•Yes — 
I<  ND 


•{Are  there  more  detectors? 


RETURN' 
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5.1.2  Subroutine  RELC0L  -  Sample  Problem  Number  1 

This  subroutine  estimates  the  iluence  from  the  density  of 

2 

collisions.  For  comparison  purposes,  4frR  times  the  fluence  is  stored 

where  R.  is  the  radius  to  each  collision  site.  Thus,  the  contribution  is 
i 


WTBC  =  particle  weight  before  collision, 

Is  =  cross  section  for  the  group  (IG0)  of  the  incoming 
particle. 


A  detector  number  from  NDC  +  1  to  ND  is  chosen  depending 
on  the  value  of  Ri  . 

Called  from:  BANKR(5) 

Subroutines  called:  FLUXST,  NSIGTA 

Commons  required:  PDET,  NUTR0N,  USER,  DETCUT 

Variables  required: 

X,  Y,  Z,  IG0,  NMED,  AGE,  WTBC  from  NUTR0N  common 
TS1G  and  detector  radii. 

Significant  internal  variables:- 

R  -  radial  distance  to  collision  site, 

C0N  -  fluence  estimate. 
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Subroutine  RELC0L 


\ 

l 
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5. 1.  3  Subroutine  BDRYX  -  Sample  Problem  Number  1  ! 

I 

This  subroutine  is  called  whenever  the  partible  in  the  random 
walk  encounters  a  change  in  geometry  media.  If  the  source-to-collision 
distance  corresponds  to  a  detector  position,  the  reciprocal  of  the  cosine 
of  the  angle  f^om  the  radius  vector  is  used  as  a  fluence  estimate. 


Called  from: 

BANKR(7) 

Subroutines  called: 

ERR0R 

Coriimons  required: 

Variables  required: 

USER, 

NUTR0N,  DET 

X,  Y,  Z, 

u,  V, 

W,  WATE  (from  common  NUTR0N) . 

Significant  internal  variables: 

R21 

- 

Radial  distance  to  boundary  crossing, 

CS3 

JB 

- 

99%  of  R21, 

R22 

- 

101%  of  R21. 

'  \  ' 

C0S 

- 

cosine  of  angle  between  particle  direction 
and  radius  vector, 

ABC0S 

- 

absolute  value  of  C0S, 

C0N 

- 

fluence  estimate. 

i 


\ 


Subroutine  BDRYX 


Calculate  R21,  distance  from  origin 
to  boundary  cross i ng . 


'  Does  R21  correspond  to  a  \ 

No  r  / Is  R21  greater  \ 

requested  detector  position?/- 

than  largest  \ 
detector  F 

Yes  |  (detector  indicated 
by  index  I) 


position  ? 


Calculate  C0S,  cosine  of] 
angle  between  particle 
direction  and  surface. 


ETURN j-  Yes  Is  the  cosine  exactly  zero?) 


>1.001 


No 


What  is  the  range  of  absolute^  <. 
value  of  C05?  _ J 


01 


ABC0S  = 
0.005 


ABC0S  =  |C0S| 


_ j 

L 

Call  FLUXST 
to  store 
contribution 

[increment 
-H  NN 

C0N  ---  WAT E /A B C0S (-*■ 

jRETUR^ 


No 
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5.  2  Sample  Problem  Number  2 

This  sample  problem  is  used  to  illustrate  the  next  event 
estimator.  The  neutron  and  secondary  gamma  ray  spectra  at  several 

4 

point  detectors  located  10  cm  from  a  polyethylene  slab  is  calculated. 
The  source  is  a  beam  of  neutrons  with  a  fission  spectrum  energy  dis¬ 
tribution.  The  use  of  the  point  detector  routines  to  estimate  the  energy 
dependent  fluence  is  illustrated  and  in  particular,  for  secondary  gamma 
rays  the  contribution  of  the  uncollided  or  first  flight  gamma  rays  are 
scored  in  subroutine  SGAM.  The  estimates  from  real  collisions  are 

made  in  subroutine  RELC0L.  Results  are  recorded  in  Table  XIX  in 

2 

units  of  particles  cm  source  neutron. 


TABLE  XIX. 

Results  of  Last  Flight  Estimator 
Sample  Problem 


Detector 

2 

Neutrons/cm 

Per 

Source  Neutron 

2 

Gammas/cm 

Per 

Source  Neutron 

1 

6.  24(-ll) 

2.  36(-ll) 

2 

6.  44(-10) 

4.  82(-ll) 

3 

1.  32(-9) 

4.  93(-ll) 

4 

1.  56(-9) 

4.  94(-ll) 

5 

6.  92(-l) 

2. 15(-3) 
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5.3 


Sample  Problem  Number  3 

This  sample  problem  serves  two  purposes  by  (1)  demonstrat¬ 
ing  the  use  of  the  PICTURE  code,  and  (2)  verifying  the  tracking  capability  • 
of  the  combinatorial  geometry  package  in  a  moderately  complex  three- 
dimensional  system.  In  the  PICTURE  program,  subroutine  G1  is  called 
from  MESH  in  the  exact  manner  that  it  is  called  from  G0MST  during  a 
M0RSE  calculation.  As  such,  PICTURE  works  by  tracking  rays  and 
thus  can  be  used  to  debug  a  combinatorial  geometry  input. 

The  sample  problem  itself  is  a  picture  of  a  tank.  The  tank 
model  itself  was  constructed  purely  as  an  illustration  of  the  combina¬ 
torial  geometry  and  is  in  no  way  accurate  or  detailed.  Notice  that  the 
use  of  both  the  0R  operator  and  the  ARB  body  is  demonstrated. 


Output  from  PICTURE  Sample  Problem 
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APPENDIX  A 

ENERGY  INDEXING  SCHEME 

For  most  problems  the  user  does  not  need  to  worry  about  the 
energy-indexing  scheme  used  in  the  M0RSE  Monte  Carlo  code.  Energies 
corresponding  to  the  group  boundaries  in  the  forward  group  structure 
are  input  and  the  code  takes  care  of  the  rest.  The  analysis  package 
SAMB0  performs  the  proper  bookkeeping  for  labeling.  A  diagram  for 
the  energy  group  structure  is  given  in  Fig.  Al.  However,  some  clarifi¬ 
cation  of  the  indexing  scheme  for  several  types  of  problems  is  advanta¬ 
geous.  This  clarification  can  perhaps  best  be  made  through  the  use  of 
examples.  Consider  a  five-group  coupled  neutron-gamma-ray  problem 
with  the  following  group  boundaries: 

Group  Energy  (eV)  with  NMGP  =  3,  NMTG  =  5, 

1  15 (+6)  NGPQTN  =  3,  NGPQTG  =  2 

2  5(+6) 

3  l(+6) 

EB0TN  3  (+3) 

4  10  (+6) 

5  2  (+6) 

EB0TG  4  (+5) 

There  are  several  options  which  might  be  considered;  namely, 
(1)  a  forward  coupled  problem,  (2)  a  forward  neutron-only  problem,  (3) 
a  forward  gamma -ray-only  problem,  (4)  an  adjoint  coupled  problem,  (5) 
an  adjoint  neutron-only  problem,  and  (6)  an  adjoint  gamma-ray-only 
problem.  Table  AI  gives  the  energies  as  they  are  indexed  during  a 
M0RSE  run.  Table  All  gives  the  values  of  variables  NQTl,  NQT2,  and 
NQT3  as  they  appear  in  USER  common  for  the  same  six  problems. 

Preceding  page  blank 
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NGPQT2 
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1 

Fig.  Al.  Diagram  of  energy  group  structure. 
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NGPQT3  =  NGPQT1  for  neutron  only  or  gamma  ray  only  problem. 
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APPENDIX  B 


INPUT  INSTRUCTIONS 


B.  1  M0RSE  Input  Instructions 

The  input  read 1  by  subroutine  INPUT  is  as  follows: 


CARD  A  (20A4) 

Title  card.  (Any  character  other  than  a  blank  or  alphameric  in 
column  one  will  terminate  the  job. ) 


CARD  B  (1315) 
NSTRT 
NM0ST 


NITS 

NQUIT 

NGPQTN 

NGPQTG 

NMGP* 


number  of  particles  per  batch, 

maximum  number  of  particles  allowed  for  in  the  \ 
bank(s);  may  equal  NSTRT  if  there  is  no  splitting, 
fission,  and  secondary  generation  during  ex¬ 
ecution.  If  bank  size  is  exceeded  by  more  than 
50  due  to  fission  or  secondary  gamma  ray  gen¬ 
eration  the  job  is  terminated, 

number  of  batches, 

number  of  sets  of  NITS  batches  to  be  run  without 
calling  subroutine  INPUT, 

number  of  neutron  groups  being  analyzed, 

number  of  gamma -ray  groups  being  analyzed, 

number  of  primary  particle  groups  for  which 
cross  sections  are  stored;  should  be  same  as 
NGP  (or  the  same  as  NGG  when  NGP  =  0)  on 
Card  XB  read  by  subroutine  XSEC, 


Input  and  output  logical  unit  members  are  defined  in  MAIN  and  ii  Gl. 
Unit  16  is  a  temporary  storage  unit  used  by  J0M3N. 

See  Table  BI  for  sample  input. 
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NMTG* 


NC0LTP  t 


IADJM 
MAXTiM  - 

MEDIA 

MEDALB  - 

\ 

I 


total  number  of  groups  for  which  cross  sections 
are  stored;  should  be  same  $s  NGP  +  NGG  as 
read  on  Card  XB  read  by  subroutine  XSEC, 

set  greater  than  zero  if  a  collision  tape  is  de¬ 
sired;  the  collision  tape-  is  written  by  the  user 
routine  BANKR, 

set  greater  than  zero  for  an  adjoint  problem, 

maximum  clock  time  in  minutes  allowed  for  the 
I  problem  to  be  on  the  computer, 

number  of  cross-section  media;  should  agree  | 
w’t h  NMED  cn  Card  XB  read  by  subroutine  XSEC, 

albedo  scattering  medium  is  absolute  value  oi 
MEDALB;  if 

=  0,  no  albedo  information  to  be  read  in, 

<0,  albedo  only  problem— no  cross  sections 
are  to  be  read, 

>  0,  coupled  albedo  and  transport  problem. 


CARD  C  (415,  5E10.  5) 


IS0UR 

NGPFS 

ISBIAS 


source  energy  group  if  >  0,  if  IS0UR  £  0, 

S0RIN  is  called  for  input  of  Cards  El  and  E2, 

rumlber  of  groups  for  which  the  source  spectrum 
is  to  be  defined.  If  IS0UR  £  0,  NGPFS  *  2, 

-  \  no  source  energy  biasing  if  set  equal  to  zero; 
otherwise  the  source  energy  is  to  be  biased, 
and  Cards  E2  are  required, 


N0TUSD 

WTSTRT 


an  unused  variable, 

! 

weight  assigned  to  each  source  particle, 


EB0TN 


lower  energy  limit  of  lowest  neutron  group  (eV) 
(group  NMGP), 


EB0TG 


lower  energy  limit  of  lowest  gamma-ray  group 
(eV)  (group  NMTG),  ; 


\ 


* 

See  Table  BI  for  sample  input; 

\ 
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I  *  W  3 VBfZ  »  f* 


^1111  ijji|i)fjp,^»^.i|i4vij^^w^^wu#-w  uupi  m 'fl 


\ 


TCUT  -  age  in  sec  at  which  particles  are  retired;  if 
TCUT  =  0,  no  time  kill  is  performed, 

VELTH  -  velocity  of  group  NMGP  when  NGPQTN  >  0;  i.  e. , 
thermal* -neutron  velocity  (cm/sec). 


CARD  D  (7E10. 4) 


XSTRT 

YSTRT 

ZSTRT 


coordinates  for  soured  particles, 


^GSTRT 

UINP 

VINP 

WINP 


starting  age  for  source  particles, 

source  particle  direction  cosines ;  if  all  are 
zero,  isotropic  directions  are  chosen. 


\ 

Source  data  on  Cards  C  and  D,  may  be  overridden  by  changes 
in  subroutine  S0URCE1. 


\ 


CARDS  El  (7E10. 4)  (omit  if  IS0UR  on  Card  C  >  0) 

NGPFS  values  of  FS,  where  FS  equal  the  unnormalized  fraction 
of  source  particles  in  each  group. 


CARDS?  E2  (7E10.  4)  (omit  if  IS0UR  >  0  or  if  IS0UR  *  0  and  LSBIAS  =  0) 

If  ISBIAS  >0,  NGPFS  values  of  BFS,  the  relative  importance  of 
a  source  in  group  I,  are  required. 

\  CARDS  F  (7E10.  4) 

NMTG  values  of  ENER,  the  energies  (in  eV)  at  the  upper  edge  of 
the  energy  group  boundaries.  (Note:  The  lower  energies  of 
groups  NMGP  and  NMTG  were  read  on  Card  C. ) 


CARD  G  (215,  5X,  3611,  1311)  (omit  if  NC0LTP  on  Card  B  £  0) 

NHISTR  -  logical  tape  number  for  the  first  collision  tape, 

NHISMX  -  the  highest  logical  number  that  a  collision  tape 

I  may  be  assigned, 
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NBIND(J),  J-l,  36  -  an  index  to  indicate  the  collision  para¬ 

meters  to  be  written  on  tape, 

NC0LLS(J),  J=l,  13  -  an  index  to  indicate  the  types  of  collisions 

to  be  put  on  tape.  (See  Table  BII  for  in¬ 
formation  concerning  NBIND's  and 
NC0LLS. ) 

CARD  H  (Z12) 

RAND0M  -  starting  random  number. 

CARD  I  (715) 

NSPLT  -  index  indicating  that  splitting  is  allowed  if  >  0, 

NKILL  -  index  indicating  that  Russian  roulette  is  allowed 

if  >0, 

NPAST  -  index  indicating  that  exponential  transform  is  in¬ 
voked  if  >  0  (subroutine  DIREC  required), 

N0LEAK  -  index  indicating  that  non-leakage  is  invoked  if 

>0, 

IEBIAS  -  index  indicating  that  energy  biasing  is  allowed 
if  >0, 

MXREG^  -  number  of  regions  described  by  geometry  input 
(’-i.11  be  set  to  one  if  *  0), 

MAXGP  -  group  number  of  last  group  for  which  Russian 
roulette  or  splitting  or  exponential  transform  is 
to  be  performed. 


t 


I'  ENDRUN  is  used,  a  data  array  relating  media  numbers  to  region 
numbers  must  be  given  in  a  data  statement  in  ENDRUN. 


li'jjM'JM 


CARD  J  (615,  4E10.  5)  (omit  if  NSPLT  +  NKILL  +  NPAST  =  0) 

from  energy  group  NGPl  to  energy  group  NGP2, 
inclusive,  in  steps  of  NDG  and  from  region  NRG1 
to  NRG2,  inclusive,  in  steps  of  NDRG,  the  follow¬ 
ing  weight  standards  and  path -stretching  para¬ 
meters  are  assigned.  If  NGPl  =  0,  groups  1  to 
MAXGP  will  be  used;  if  NRG1  =  0,  regions  1  to 
MXREG  will  be  used  (both  in  steps  of  one).  Us¬ 
ually  NDG  =  1  and  NDRG  =  1. 

weight  above  which  splitting  will  occur, 

weight  below  which  Russian  roulette  is  played, 

weight  given  those  particles  surviving  Russian 
roulette, 

path-length  stretching  par .  meters  for  use  in  ex  ¬ 
ponential  transform  (usually  0  £  PATH  <  1). 

The  above  information  is  repeated  until  data  for  all  groups  and 
regions  are  input. 

End  Cards  J  with  negative  value  of  NGPl  (ex. ,  -1  in  columns  4 
and  5). 

CARDS  K  (7E10. 4)  (omit  if  IEBIAS  on  Card  I  *  0) 

((EPR0B(IG,  NREG),  IG  =  1,  NMTG),  NREG  =  1,  MXREG) 

Values  of  the  relative  energy  importance  of  particles  leaving  a 
collision  in  region  NREG.  Input  for  each  region  must  start  on  a 
new  card. 

CARD  L  (415) 

NS0UR  -  set  £  0  for  a  fixed  source  problem;  otherwise  the 
source  is  from  fissions  generated  in  a  previous 
batch, 

MFISTP  -  index  for  fission  problem,  if  £  0  no  fissions  are 
allowed, 

NKCALC  -  the  number  of  the  first  batch  to  be  included  in  the 
estimate  of  k;  if  *  0  no  estimate  of  k  is  made, 


NGPl 

NDG 

NGP2 

NRG1 

NDRG 

NRG2 


WTHIH 

WTL0W 

WTAVE 


PATH 
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N0RMF 


the  weight  standards  and  fission  weights  are  un¬ 
changed  if  £  0;  otherwise  fission  weights  will  be 
multiplied,  at  the  end  of  each  batch,  by  the  latest 
estimate  of  k  and  the  weight  standards  are  mul¬ 
tiplied  by  the  ratio  of  fission  weights  produced  in 
previous  batch  to  the  average  starting  weight  for 
the  previous  batch.  For  time -dependent  decay¬ 
ing  systems,  N0RMF  should  be  >  0. 


CARDS  M  (7E10. 4)  (omit  if  MFISTP  on  Card  L  *  0) 

(FWL0(I),  1=1,  MXREG)  values  of  the  weight  to  be  assigned  to 
fission  neutrons. 


CARDS  N  (7E10. 4)  (omit  if  MFISTP  on  Card  L^O) 

(FSE(IG,  IMED),  IG  =  1,  NMGP),  IMED  =  1,  MEDIA)  the  fraction 
of  fission -induced  source  particles  in  group  IG  and  medium  IMED. 

Note:  Input  for  each  medium  must  start  on  a  new  card. 


CARDS  0  (7E10.  5)  (omit  if  NGPQTN  =  0  or  NGPQTG  =  0,  i.  e. ,  include 
if  coupled  neutron-gamma-ray  problem) 

((GWL0(IG,  NREG),  IG  =  1,  NMGP  or  NMTG  -  NMGP),  NREG  =  1, 


MXREG) 


average  number  of  secondary  particles  to  be 
generated.  NMGP  groups  are  read  for  each 
region  in  a  forward  problem  and  NMVG  - 
NMGP  for  an  adjoint.  Input  for  each  region 
must  start  on  a  new  card. 
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TABLE  BI. 


Sample  Group  Input  Numbers  for  Some 
Representative  Problems* 


Case  A  -  Neutron  Only  Cross  Sections  (22  groups) 

Case  B  -  Gamma -Ray  Onlv  Cross  Sections  (18  groups) 

Case  C  -  Neutron-Gamma -Ray -Coupled  Cross  Sections  (22-18  groups) 


Input  Variable 

Problem  Type 

Case  A 

Top  14  groups 

Case  B 

Top  17  groups 

Case  C 

Neutrons  Only 

Top  14  groups 

Case  C 

Gamma  Rays  Only 
Top  17  groups 

Case  C 

Neutron-Gamma 
Top  14  Neutron, 
Top  17  Gamma 
groups 

M0RSE  Input: 

NGPQTN 

14 

0 

14 

0 

14 

NGPQTG 

0 

17 

0 

17 

17 

CARD  B, 

NMGP 

22 

18 

22 

18 

22 

Variables 

NMTG 

22 

18 

22 

18 

40 

NGP 

22 

18 

22 

0 

22 

NGG 

0 

0 

0 

18 

18 

CARD  XB 

INGP 

22 

18 

40 

40 

40 

Variables 

T 


For  cross  sections  with  full  downscatter,  NDS  =  NGP,  NDSG  =  NGG, 
INDS  =  iNGP,  and  ITBL  =  number  of  downscatters  +  number  of  up- 
scatters  +  3.  Usually,  ISGG  =  number  of  upscatters  +  4. 
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TABLE  BII. 


BANKR  Arguments  (NC0LLS) 


BANKR 

Argument 

Called 

from 

Location  of  Call  in  Walk 

-1 

M0RSE 

After  call  to  INPUT— to  set  parameters 
for  new  problem 

-2 

M0RSE 

At  the  beginning  of  each  batch  of  NSTRT 
particles 

-3 

M0RSE 

At  the  end  of  each  batch  of  NSTRT 
particles 

-4 

M0RSE 

At  the  end  of  each  set  of  NITS  batches— 
a  new  problem  is  about  to  begin 

1 

MS0UR 

After  a  source  event 

2 

TESTW 

After  a  splitting  has  occurred— com¬ 
mented  in  column  1 

3 

FPR0B 

After  a  fission  has  occurred 

4 

GST0RE 

After  a  secondary  particle  has  been 
generated 

5 

M0RSE 

After  a  real  collision  has  occurred— 
post -collision  parameters  are  avail¬ 
able 

6 

M0RSE 

After  an  albedo  collision  has  occurred — 
post -collision  parameters  are  available 

7 

NXTC0L 

After  a  boundary  crossing  occurs  (the 
track  has  encountered  a  new  geometry 
medium  other  than  the  albedo  or  void 
media) 

8 

NXTC0L 

After  an  escape  occurs  (the  geometry 
has  encountered  medium  zero) 

9 

M0RSE 

After  the  post-collision  energy  group 

exceeds  the  maximum  desired— com¬ 
mented  in  column  1 


r-i'T'V 


TABLE  BH.  (Cont’d.) 


BANKR 

Argument 

Called 

from 

Location  of  Call  in  Walk 

10 

M0RSE 

After  the  maximum  chronological  age 
has  been  exceeded— commented  in 
column  1 

11 

TESTW 

After  a  Russian  roulette  kill  occurs— 
commented  in  column  1 

12 

TESTW 

After  a  Russian  roulette  survival 
occurs— commented  in  column  1 

13 

GST0RE 

After  a  secondary  particle  has  been 
generated  but  no  room  in  the  bank  is 
available— commented  in  column  1. 
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B.  2  05R  Geometry  Input  Instructions 

Read  by  subroutine  J0MIN  and  the  specific  input  depends  on  the 
geometry  packages  used. 


B.  2. 1  Spherical  GE0M 
CARD  GA  (15,  DIO.  5) 

MED  -  medium  number  interior  to  R(>  0), 

R  -  outer  radius  of  sphere  or  spherical  shell  contain¬ 

ing  MED. 

Repeat  Card  GA  for  all  radii  (s  20)  in  increasing  order.  End 
Card  GA  input  with  blank  card. 


CARD  GB  (DIO.  5) 

R  -  region  radius  of  sphere  or  spherical  shell  con¬ 

taining  regions.  Region  numbers  are  assigned 
in  consecutive  order  starting  with  1,  and  R  must 
be  in  increasing  order. 

Repeat  Card  GB  for  all  radii  (s  20)  in  increasing  order.  End  Card 
GB  input  with  blank  card.  If  no  regions  are  desired,  a  blank  card 
must  be  used  to  signal  no  region  geometry. 

B.2.2  Slab  GE0M 

CARD  GA  (15,  DIO.  5) 

MED  -  medium  with  Z  as  lower  bound  (>  0), 

Z  -  lower  limit  of  medium  MED. 

Repeat  Card  GA  for  all  boundaries  with  the  last  card  containing 
MED  =  0  and  the  boundary  of  the  system. 
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CARD  GB  (DIO.  5) 

Z  -  lower  limit  of  region  boundary.  Region  numbers 

are  assigned  in  consecutive  order  starting  with 
1,  and  Z  must  be  in  increasing  order. 

Repeat  Card  GB  for  all  region  boundaries. 

End  Card  GB  input  with  a  blank  card.  If  no  region  geometry  is 
desired,  a  blank  card  is  required. 


CARD  GC  (4D10.  5) 

XL 

lower 

XU 

upper 

YL 

lower 

YU 

upper 

boundary  of  system  in  X  direction, 
boundary  of  system  in  X  direction, 
boundary  of  system  in  Y  direction, 
boundary  of  system  in  Y  direction. 


B.  2. 3  Cylindrical  GE0M 
CARD  GA  (15,  5X,  A8) 

NREGIN  -  flag  to  indicate  material  media  (=2)  or  both 
region  and  material  media  (=1), 

SEX  -  sex  of  programmer. 


CARD  GB  (E10.  5) 

R  -  radii  of  the  cylindrical  shells  describing  the 

material  media  in  ascending  order. 

Repeat  Card  GB  until  all  radii  have  been  input. 

End  Card  GB  input  with  a  blank  card. 
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CARD  GC  (ElO.  5,  1215/815) 

H  -  upper  height  of  medium  M(I)  (>0),  cylinders  as¬ 

sumed  to  start  at  H  =  0, 

M(I)  -  media  for  the  cylindrical  shells  for  this  height. 

Repeat  Card  GC  until  all  height  intervals  have  been  input. 

End  Card  GC  input  with  a  blank  card  or  if  there  are  more  than  12 
radial  intervals,  2  blank  cards. 


CARD  GC  (E10.  5)  (omit  if  NREGIN  =  2) 

RG  -  radii  of  the  cylindrical  shells  describing  the 

region  geometry  in  ascending  order. 

Repeat  Card  GD  until  all  region  geometry  has  been  input. 

End  Card  GD  with  a  blank  card. 

CARD  GE  (E10.  5,  1215/815)  (omit  if  NREGIN  =  2) 

HG  -  upper  height  of  region  MG (I), 

MG(I)  -  region  numbers  for  the  cylindrical  shells  for  this 
height. 

Repeat  Card  GE  until  all  height  intervals  have  been  input. 

End  Card  GE  input  with  a  blank  card  or  if  there  are  more  than  12 
radial  intervals,  2  blank  cards. 


B.  2.4  General  GE0M 


CARD  GA  (15,  5X,  A6, 
NSTAT 

SEX 

STATUS 


IX,  A7)  hollerith  left  adjusted 

flag  to  indicate  material  media  only  if  2  and  both 
region  and  material  media  if  1, 

sex  of  the  programmer  (select  one  from  MALE, 
FEMALE,  or  blank  indicating  uncertain), 

marital  status  of  programmer. 
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CARD  GB  (2A4,  A3,  5(D10.  5,  Al)) 

DUMMY(3)  -  bollerith  characters  not  used, 

FIN(I)  -  zone  boundaries  increasing  order  along  the  x 
axis, 

BCD(I)  -  flag  to  indicate  end  of  input  if  blank,  comma 
means  to  continue. 

Repeat  in  format  (6(D10.  5,  Al))if  more  than  five  boundaries 
along  the  x  axis  are  needed. 


CARD  GC  -  same  as  CARD  GB  except  for  y  axis. 

Repeat  in  format  (6(D10.  5,  Al))  if  more  than  five  boundaries 
along  the  y  axis  are  needed. 


CARD  GD  -  same  as  Card  GB  except  for  z  axis. 

Repeat  in  format  (6(D10.  5,  Al))  if  more  than  five  boundaries 
along  the  z  axis  are  needed. 


CARD  GE  (A4,  A2,  315) 

BCD1  -  hollerith  Z0NE, 

BCD2  -  dummy, 

integers  which  specify  the  zone  as  being  the 
NXZN0th  zone  in  the  x  direction,  NYZN0th  zone 
in  the  y  direction,  and  NZZN0th  zone  in  the  z 
direction. 


CARD  GF  (2A4,  A3,  5(D10.  5,  Al)) 

DUMMY(3)  -  hollerith  characters  not  used. 

FIN (I)  -  block  boundaries  in  increasing  order  along  the 

x  axis, 


NXZN0 

NYXN0 

NZZN0 
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BCD(I)  -  flag  to  indicate  end  of  input  if  blank,  comma 
means  to  continue. 

Repeat  in  format  (6(D10.  5,  Al))  if  more  than  five  boundaries 
along  the  x  axis  are  needed. 

CARD  GG  -  same  as  Card  GF  except  for  y  axis. 

Repeat  in  format  (6(D10.  5,  Al))  if  more  than  five  boundaries 
along  the  y  axis  are  needed. 

CARD  GH  -  same  as  Card  GF  except  for  z  axis. 

Repeat  in  7ormat  (6 (DIO.  5,  Al))  if  more  than  five  boundaries 
along  the  z  axis  are  needed. 

CARDS  GI  to  G0  describe  the  geometry  for  a  block  and  must  be  included 
for  each  block  in  me  zone. 

CARD  GI  (A4,  A2,  315) 

BCD1  -  hollerith  BL0C, 

BCD2  -  dummy, 

NXBND  )  integers  which  specify  the  block  as  being  the 
NYBND  |  -  NXBNDth  in  the  x  direction,  the  NYBNDfh  in  the 
NZBND  )  y  direction,  and  the  NZBND  in  the  z  direction. 

CARD  GJ  (3A4,  10(15,  Al)) 

NAM2  -  hollerith  MEDI, 

DUM{2)  -  dummy, 

INP(I)  -  a  list  of  media  sector  by  sector  in  the  block, 

BCD(I)  -  flag  to  indicate  end  of  input  if  blank,  a  comma 

means  to  continue. 

Continuation  with  12(15,  Al). 
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CARD  GK  (3A4,  10(15,  Al)) 

NAM2  -  hollerith  SURF, 

DUM(2)  -  dummy, 

INP(I)  -  a  list  of  quadratic  surfaces  appearing  in  the 
block  Numbers  must  appear  in  the  order  the 
surfaces  are  described  on  Card  GQ, 

SCD(I)  -  flag  to  indicate  end  of  input  if  blank,  a  comma 
means  to  continue. 

Continuation  of  Card  GK  in  12(15,  Al)  format  is  permissible. 

CARD  GL  (A4,  A2,  1813) 

SI  -  hollerith  SECT, 

DUM  -  dummy, 

IND(I)  -  the  designation  of  each  sector  wliich  describes 

the  position  of  the  sector  relative  to  quadratic 
surfaces, 

+1:  sector  is  on  positive  side  of  surface, 

-1:  sector  is  on  negative  side  of  surface, 

0:  surface  is  not  needed  to  define  sector. 

There  must  be  a  Card  GL  for  each  sector  and  references  to 
quadratic  surfaces  must  be  in  same  order  as  they  are  listed 
on  Card  GQ. 

CARD  GM  (3A4,  10(15,  Al))  (omit  if  NSTAT  on  Card  GA  =  2) 

NAM2  -  hollerith  REGI, 

DUM(2)  -  dummy, 

INPOO  -  a  list  of  regions  sector  by  sector  in  the  block, 

BCD(I)  -  flag  to  indicate  end  of  input  if  blank,  a  comma 
means  to  continue. 

Continuation  with  12(15,  Al)  format  is  permissible. 
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CARD  GN  (3A4,  10(15,  Al))  (omit  if  NSTAT  on  Card  GA  =  2) 


NAM2 

DUM(2)  , 

INP(I) 

BCD(I) 


hollerith  SURF, 

i  \ 

dummy, 

same  as  for  Card  GK  except  for  region  input  in¬ 
stead  of  material  inpu^. 


CARD  GO  (A4,  A2,  1813)  (omit  if  NSTAT  on  Card  GA  =  2) 


DUM 


hollerith  SECT, 
dummy,  1 

same  as  for  Card  GL  except  for  region  input  in- 


IND(I)  -  same  as  for  Card  GL  e: 

stead  of  material  input. 

Repeat\Cards  GI  to  G0  for  each  block. 


CARD  GP  (15,  16A4,  A2)  ; 


N0BD 


DljM(I) 


total  number  of  quadratic  surfaces  in  the  entire 
system,  '■ 

hollerith  characters  ignored  by  the  code.  (Help¬ 
ful  in  identifying  input  at  a  later  time. )  1 


CARD  GQ  (4  (DIO.  5j  A4,  IX,  Al)) 


C0F(J) 

BCD1(J) 

i 

i 

BCD2(J) 


coefficient  of  the  term, 

\ 

hollerith  indicating  which  term  of  the  equation. 
XJQ,  YSQ,  ZSQ,  XZ,  YX,  YZ,  XY,  ZX,  YZ, 
X,  Y,  Z,  or  blank  are  the  possibilities, . 


BCD2(J)  -  a  flag  which  indicates  the  quadratic  equation 
continues.  Any  non-blank  character  ends  fhe 

f|eld.  The  next  function  .must  start  on  row  card, 

\ 

Repeat  Caras  GQ  until  all  surfaces  have  been  described. 

A  sample  of  the  input  is  shown  in  Table  BUI. 


144 


ujwi'wwm1,1 '  «  1 1  y  *r 


B.3  Combinatorial  Geometry  Input  Instructions 

The  combinatorial  geometry  input  data  is  read  by  the  J0MIN 
subroutine,  except  for  the  region  volumes  VN0R(I),  which  are  read  by 
the  GTVLIN  subroutine  whenever  IV0PT  =  3.  For  clarity  of  terminology, 
the  terms  ’’regions”  and  ’’media"  have  essentially  the  same  meaning  as 
in  the  05R  Geometry  Package,  but  are  constructed  in  a  different  manner. 
The  term  "zone”  is  the  same  as  the  "region”  as  defined  in  the  original 
combinatorial  geometry  package.  The  term  "body"  has  the  same  mean¬ 
ing  as  in  the  original  combinatorial  geometry  package. 


CARD  CGA  (215,  10X,  10A6) 

IV0PT  -  option  which  defines  the  method  by  which  region 
volumes  are  determined;  if 

IV0PT  =  0,  volumes  set  equal  to  1. , 

IV0PT  =  1,  concentric  sphere  volumes  are  cal¬ 
culated 

IV0PT  =  2,  slab  volumes  (1-dim. )  are  calculated, 
IV0PT  =  3,  volumes  are  input  by  card, 


IDBG  -  if  IDBG  >0,  subroutine  PR  is  called  to  print  re¬ 
sults  of  combinatorial  geometry  calculations 
during  execution.  Use  only  for  debugging, 

JTY  -  alphanumeric  title  for  geometry  input  (columns 

21-80). 


CARDS  CGB  (2X,  A3,  IX,  14,  6E10.  3) 

One  set  of  CGB  cards  is  required  for  each  body  and  for  the  END 
card  (see  Table  B  IV).  Leave  columns  1-6  blank  on  all  continu¬ 
ation  cards. 


* 


Not  operational. 


ITYPE 

IALP 


FPD(I) 


specifies  body  type  or  END  to  terminate  reading 
of  body  data  (for  example  B0X,  RPP,  ARB,  etc.). 
Leave  blank  for  continuation  cards, 

body  number  assigned  by  user  (all  input  body 
numbers  must  form  a  sequence  set  beginning  at 
1).  If  left  blank,  numbers  are  assigned  sequen¬ 
tially.  Either  assign  all  or  none  of  the  numbers. 
Leave  blank  for  continuation  cards, 

real  data  required  for  the  given  body  as  shown 
in  Table  B  IV. 


CARDS  CGC  (2X,  A3,  15,  9(A2,  15)) 

Input  zone  specification  cards.  One  set  of  cards  required  for 

each  input  zone,  with  input  zone  numbers  being  assigned  sequen¬ 
tially. 

IALP  -  IALP  must  be  a  nonblank  for  the  first  card  of 
each  set  of  cards  defining  an  input  zone.  If 
IALP  is  blank,  this  card  is  treated  as  a  con¬ 
tinuation  of  the  previous  zone  card, 

IALP  =  END  denotes  the  end  of  zone  description. 

NAZ  -  total  number  of  zones  that  can  be  entered  upon 

leaving  any  of  the  bodies  defined  for  this  input 
region  (some  zones  may  be  counted  more  than 
once).  Leave  blank  for  continuation  cards  for  a 
given  zone.  (If  NAZ  £  Oon  the  first  card  of  the 
zone  card  set,  then  it  is  set  to  5).  This  is  used 
to  allocate  blank  common, 

Alternate  nBIAS(I)  and  JTY(I)  for  all  bodies  de¬ 
fining  this  input  zone. 

IIBIAS(I)  -  specify  the  "OR”  operator  if  required  for  the 
JTY(I)  body, 

JTY(I)  -  body  number  with  the  (+)  or  (-)  sign  as  required 
for  the  zone  description 
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CARDS  CGD  (1415) 


MRIZ(I)  -  MRIZ(I)  is  the  region  number  in  which  the  "I  h" 
input  zone  is  contained  (I  =  1,  to  the  number  of 
input  zones).  Region  numbers  must  be  sequen¬ 
tially  defined  from  1. 


CARDS  CGE  (1415) 

MMIZ(I)  -  MMIZ(I)  is  the  medium  number  in  which  the 
tqthtf  input  zone  is  contained  (I  =  1,  to  the 
number  of  input  zones).  Medium  numbers  must 
be  sequentially  defined  from  1. 


CARDS  CGF  (7E10.  5)  (omit  if  IV0PT  4  3) 

VN0R(I)  -  volume  of  the  "Ith"  region  (I  =  1  to  MXREG,  the 
number  of  regions). 

Note:  If  ENDRUN  is  used  to  obtain  collision  density  and  track 
length  per  unit  volume  estimate  of  fluence,  then  a  data 
statement  in  ENDRUN  must  give  a  relationship  between 
region  and  media.  In  this  case  only  one  medium  may  be 
in  a  region. 
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Card  Columns  ITYPE  IALP  Real  Data  Defining  Particular  Body  Number  of 

Body  Type  3-5  7-10  11-20  21-30  31-40  41-50  51-60  61-70  Cards  Needed 
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B.  4  M0RSEC— Cross -Section  Module  Input  Instructions 


CARD  XA  (20A4) 

Title  card  for  cross  sections.  This  title  is  also  written  on  tape 
if  a  processed  tape  is  written;  therefore,  it  is  suggested  that  the 
title  be  definitive. 


CARD  XB  (1315) 


Column 

Name 

Description 

5 

NGP 

the  number  of  primary  groups  for  which 
there  are  cross  sections  to  be  stored. 

Should  be  same  as  NMGP  input  in  M0RSE, 

10 

NDS 

number  of  primary  downscatters  for  NGP 
(usually  NGP), 

15 

NGG 

number  of  secondary  groups  for  which  there 
are  cross  sections  to  be  stored, 

20 

NDSG 

number  of  secondary  downscatters  for  NGG 
(usually  NGG), 

25 

INGP 

total  number  of  groups  for  which  cross 
sections  are  to  be  input. 

30 

ITBL 

table  length,  i.  e. ,  the  number  of  cross 
sections  for  each  group  (usually  equal  to 
number  of  downscatters  +  number  of  up- 
scatters  +  3), 

35 

ISGG 

location  of  within -group  scattering  cross 
sections  (usually  equal  to  number  of  up- 
scatters  +  4), 

40 

NMED 

number  of  media  for  which  cross  sections 
are  to  be  stored- -should  be  same  as  MEDIA 
input  in  M0RSE, 

45 

NELEM 

number  of  elements  for  which  cross  sections 
are  to  be  read, 

CARD  XB  (Cont’d. ) 
Column  Name 


Description 


50 

NMDC 

number  of  mixing  operations  (elements 
times  density  operations)  to  be  performed 
(must  be  ^  1), 

55 

NC0EF 

number  of  coefficients  for  each  element, 
including  Pq, 

60 

NSCT 

number  of  discrete  angles  (usually 

NC0EF/Wal), 

65 

ISTAT 

flag  to  store  Legendre  coefficients  if 
greater  than  zero. 

CARD  XC  (1115) 


Column 

Name 

Description 

5 

IRDSG  + 

switch  to  print  the  cross  sections  as  they 
are  read  if  >  0,  if  <  0  card  sequence  is  not 
checked, 

10 

ISTR+ 

switch  to  print  cross  sections  as  they  are 
stored  if  >  0, 

15 

IFMU+ 

switch  to  print  intermediate  results  of  p’s 
calculation  if  >  0, 

20 

IM0Mf  ' 

switch  to  print  moments  of  angular  dis¬ 
tribution  if  >  0, 

25 

IPRINf 

switch  to  print  angles  and  probabilities  if 
>0, 

30 

IPUN+ 

switch  to  print  results  of  bad  Legendre  co¬ 
efficients  if  >  0, 

35 

IDTFf 

switch  to  signal  that  input  format  is  DTF-IV 
format  if  >  0;  otherwise,  ANISN  format  is 
assumed, 

CARD  XC  (Cont’d. ) 


Column  Name 


Description 


40  IXTAPE  -  logical  tape  unit  if  binary  cross-section 

tape,  set  equal  to  0  if  cross  sections  are 
from  cards.  If  negative,  then  the  pro¬ 
cessed  cross  sections  and  ether  necessary 
data  from  a  previous  run  will  be  read;  in 
this  case  (IXTAPE  <  0)  no  cross  sections 
from  cards  and  no  mixing  cards  may  be  in¬ 
put.  The  absolute  value  of  IXTAPE  is  the 
logical  tape  unit, 

45  JXTAPE  -  logical  tape  unit  of  a  processed  cross- 

section  tape  to  be  written.  This  processed 
tape  will  contain  the  title  card,  the  vari¬ 
ables  from  common  L0CSIG  and  the  per¬ 
tinent  cross  sections  from  blank  common, 

50  I06RT  -  logical  tape  unit  of  a  point  cross-section 

tape  in  06R  format, 

55  IGQPT  -  last  group  (M0RSE  multigroup  structure) 

for  which  the  06R  point  cross  sections  are 
to  be  used  NMGP). 


CARD  XD  (1415)  (omit  if  IXTAPE  £  0) 

Element  identifiers  for  cross -section- tape.  If  element  identifiers 
are  in  same  order  as  elements  on  tape,  the  efficiency  of  the  code 
is  increased  due  to  fewer  tape  rewinds. 


CARD  XE  (omit  if  IXTAPE  i  0) 

Cross  sections  in  ANISN  format  if  IDTF  *  0,  otherwise,  DTF-IV 
format.  Cross  sections  for  INGP  groups  with  a  table  length  ITBL 
for  NELEM  elements  each  with  NC0EF  coefficients. 
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CARDS  XF  (215,  E10.  5)  (NMIX  cards  are  required,  omit  if 

*  IXTAPE  <0) 


KM  -  medium  number  (media  numbers  from  1  to 

MEDIA,  see  Card  B,  must  appear  on  some  XF 
card), 

KE  -  element  number  occurring  in  medium  KM  (neg¬ 

ative  value  indicates  last  mixing  operation  for 
that  medium  and  at  least  one  negative  value  is 
required  for  each  medium), 

RH0  -  density  of  element  KE  in  medium  KM  in  units  of 

atoms/(barn  cm). 


CARDS  XG  (15)  (omit  if  I06RT  *  0) 

NXPM  -  number  of  point  cross-section  sets  per  medium 
found  on  an  06R  tape, 

=  1,  total  cross  section  only, 

=  2,  total  and  scattering  cross  section, 

=  3,  total,  scattering,  and  v  *fission  cross 
section. 

XCHEKR  Card  (415)  Cross  sections  and  cross-section  input  data  may 
be  checked  independently  of  M0RSE  utilizing  XCHEKR.  (5)  The 
input  to  XCHEKR  consists  of  the  cross-section  cards  XA  through 
XG  preceded  by  a  card  as  follows: 

IADJM  -  set  greater  than  zero  for  an  adjoint  problem, 

MEDIA  -  number  of  cross-section  media;  should  equal 

NMED  on  Card  XB, 

NMGP  -  number  of  primary  particle  energy  groups  for 
which  cross  sections  are  to  be  stored;  should 
equal  NGP  (or  NGG  if  NGP  =  0)  on  Card  XB, 

NMTG  -  total  number  of  energy  groups  for  which  cross 
section  are  to  be  stored.  Should  be  equal  to 
NGP  +  NGG  on  Card  XB. 
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B.  5  SAMB0  Analysis  Input  Instructions 

The  following  data  are  read  from  cards  by  SC0RIN: 


CARD  AA  (20A4) 

Title  information— will  be  immediately  output. 


CARD  AB  (815) 

ND  -  number  of  detectors  (set  =1  if  £0), 

NNE  -  number  of  primary  particle  energy  bins  to  be 

used,  (must  be  ^  NE), 

NE  -  total  number  of  energy  bins  (set  =0  if  *  1), 

NT  -  number  of  time  bins  for  each  detector  (may  be 

negative,  in  which  case  I  NT  j  values  are  to  be 
read  and  used  for  every  detector)  (set  =0  if 
:  NT  I  S  !), 

NA  -  number  of  angle  bins  (set  =0  if  £  1), 

NRESP  -  number  of  energy -dependent  response  functions 

to  be  used  (set  =1  if  0), 

NEX  -  number  of  extra  arrays  of  size  NMTG  to  be  set 

aside  (useful,  for  example,  as  a  place  to  store 
an  array  of  group -to -group  transfer  probabilities 
for  estimator  routines.  If  the  subroutine 
ENDRUN  which  outputs  fluence  estimates  from 
collision  and  track  lengths  is  used,  then  this 
number  must  be  at  least  MXREG  +  2  (see  Card  I 
for  MXREG), 

NEXND  -  number  of  extra  arrays  of  size  ND  to  be  set 

aside  (useful,  for  example,  as  a  place  to  store 
detector -dependent  counters). 
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CARDS  AC  (3E10. 4)  (ND  cards  will  be  read) 

X,  Y,  Z  -  detector  location.  (If  other  than  point  detectors 
are  desired,  the  point  locations  must  still  be  in¬ 
put  and  can  be  combined  with  additional  data 
built  in  to  user  routines  to  fully  define  each  de¬ 
tector.  ) 

Note  that  the  distanc  c.  between  the  above  points  and  the  XSTRT, 
YSTRT,  ZSTRT  values  and  the  initial  age,  AGSTRT,  will  be 
used  to  define  the  lower  limit  of  the  first  time  bin. 


CARD  AD  (20A4) 

Title  or  units  for  total  responses  for  all  detectors.  Will  be  used 
in  columns  54  through  133  of  the  title  for  the  print  of  these  arrays. 


CARD  AE  (20A4) 

Title  or  units  for  each  total  response  for  all  detectors. 


CARDS  AF  (7E10. 4) 

Response  function  values.  NMTG  values  will  be  read  in  each  set 
of  AF  cards.  Input  order  is  from  energy  group  1  to  NMTG  (order 
of  decreasing  energy). 

Note:  Cards  AE  and  AF  are  read  in  the  following  order:  AEj, 
AFj, . . . ,  AFn,  AE2,  AFj,  . . . ,  AFn,  etc.  NRESP  sets  of 
AE,  AF  cards  will  be  read. 


CARD  AG  (20A4)  (omit  if  NE  s  1) 

Units  for  energy -dependent  fluence  for  all  detectors. 
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CARDS  AH  (1415)  (omit  ii  NE  ^  1)  \ 

k  \  ! 

Energy  group  ^umbers  defining  lower  limit  of  energy  bins  (in 
order  of  mercasing  group  numberJT” The  NNE  (if  >0)  entry  must 
equal  NGPQTN;  the  NE  entry  must  be  set  to  NMGP  NGPQTG 
for  a  combined  problem,  or  else  NGPQTG  or  NGPQTN. 

\ 

CARD  AI  (20A4)  (omit  if  j  NT  U  1)  1 

!  1  \ 
Units  for  time -dependent  total  responses  for  all  detectors.  ! 

CARD  AJ  (20A4)  (omit  if  |  NT  J  £  1  or  NE  ^  1)  ' 

1  1  \ 

Units  for  time  and  energy -dependent  fluence  for  all  detectors. 

1  # 

i  > 

CARDS  AK  (7E10.  4)  (omit  if  |  NT  |  £  1) 

NT  val-es  of  upper  limits  of  time  bins  for  each  detector  (in  order 

tL1™,rtasins,h™e  and  detector  number).  •  The" values  for  each 
detector  must  start  on  a  new  card.  |  NT  [  values  only  are  read 
if  NT  is  negative.  They  are  then  used  for  every  detector. 

CAR  D  AL  (20  A4)  (omit  if  NA  si) 

l  ' 

i  i 

Units  for  angle-  and  energy-dependent  fluence  for  all  detectors. 


CARD  AM  (7E10. 4)  (omit  if  NA  £  1) 


\ 


5*  Sf  i  upper  Iimits  of  “Sle  bins  (actually  cosine  bins: 
the  WA1'*  value  must  equal  one).  . 

Following  the  input  for  the  SAMB0  analysis  module,  input  cards 
for  user -written  routines  INSC0R,  S0URCE,  and  ENDKUN. 
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\ Input  to  PICTURE 


Card  PA:  Format  (15) 


NUSE: 


The  number  of  characters  to  reald  on  Card  B  to  replace 
the  standard  values  of  ATABLE.  Leave  Card  A  blank 
and  omit  Card  B  if  the  standard  ATABLE  is  desired. 
NUSE  £  50. 


Card  PB:  Format  (50A1)  (omit  if  NUSE  =  0) 

ATABLE  (I),  I  =  1,  NUSE:  The  list  of  characters  that  are  to  be 
printed  for  each  medium.  For  medium  N,  ATABLE 
(N+l)  is  printed.  \  If  N  s  47,  ATABLE  (48)  is  printed. 
The  standard  values  of  ATABLE  are: 


Medium  Number 
0  (external  void) 
1  through  9 
10  through  35 
36  through  46 

^  47  (including 
internal  voids) 


Character  Printed 

1  through  9 

A  through  Z 

various  special, 
characters 

(blank) 


GE0M  input:  Combinatorial  geometry  input. 
C^rd  PC:  Format  (212,  18A4) 


ICNT: 


IRG: 


=  0  After  this  picture,  return  to  Card  PC  for  another 
picture  with  the  same  geometry. 

=  1  After  this  picture,  return  to  Card  PA  to  read  in  a 
ne^  GE0M  input. 

=  -1  Display  the  region  geometry. 

=  0  Display  the  zone  geometry. 

=  1  Display  the  material  geometry. 


'TITLE(I),  I  =  1,  18:  72  characters  to  be  printed  as  a  title. 


159 


5-: 


Card  PD:  Format  (6E10.  5) 


X,  Y,  and  Z  coordinates  in  the  combinatorial  geometry 
of  the  upper  left  corner  of  the  picture. 


X,  Y,  and  Z  coordinates  in  the  combinatorial  geometry 
of  the  lower  right  corner  of  the  picture. 


NOTE:  Card  PD  partially  describes  the  plane  of  the  slice  by 

defining  two  points  in  the  plane  and  designates  the  top, 
bottom,  left  and  right  sides  of  the  picture. 


6. 


Card  PE:  Format  (6E10.  5) 

ux  ) 

l  Direction  numbers  proportional  to  the  direction  cosines 
UY  /  for  the  U  axis  of  the  picture.  The  U  axis  points  down  the 

„  1  printed  page  in  the  direction  the  page  moves  through  the 

UZ  /  printer. 

vx  } 

Vy  \  Direction  numbers  for  the  V  axis  of  the  picture.  The  V 
(  axis  points  to  the  right  across  the  page. 

vz  ) 


NOTE:  Card  PE  completes  the  description  of  the  plane  of  the 

slice  by  giving  a  line  in  the  plane,  also  specifies  the 
orientation  of  the  picture  on  the  output. 


7.  Card  PF:  Format  (215,  2E10.  5) 

NU:  Number  of  intervals  to  print  along  the  U  axis  (overrides 

DELU) . 

NV:  Number  of  intervals  to  print  clong  the  Y  axis  (overrides 

DELV) . 
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DELU:  Spacing  xin  GE0M  units)  of  intervals  along  the  U  axis. 
DELV:  Spacing  (in  GE0M  units)  of  intervals  along  the  V  axis. 

NOTE:  All  four  entries  are  not  required  as  input  on  Card  PF, 

see  below  for  explanation. 


Options: 

*•  XLR=yUl=ZLR=0- 

For  this  case  NU  and  NV  must  be  specified.  In  addition,  either 
DELU  or  DELV  must  be  specified.  If  the  other  is  left  blank,  the 
code  will  produce  an  undistorted  picture.  If  both  DELU  and  DELV 
are  specified  the  picture  is  likely  to  be  distorted.  Th?'  standard 
printers  give  10  characters  to  the  inch  across  a  line  but  only  6 
lines  per  inch  down  the  page.  Because  of  this  DELV  =  .  6*DELU 
is  necessary  to  produce  an  undistorted  picture. 

2-  XLR  or  YLR  0r  ZLR  ^  °* 

If  any  one  variable  on  Card  PD  is  specified,  the  code  will  calculate 
the  others  to  produce  an  undistorted  picture. 

If  both  NU  and  DELU  (or  both  NV  and  DELV)  are  specified,  DELU 
(or  DELV)  will  be  ignored. 

The  U  and  V  axes  may  have  arbitrary  orientation.  (If  they  are  not 
orthogonal,  the  resulting  picture  will  be  distorted. )  In  Option  1, 
the  first  point  will  be  at  (X,  Y,  Z)ul>  and  the  remaining  points  in 
the  directions  and  at  the  distances  specified.  In  Option  2,  the 
range  from  Xul  to  X^r  is  -’-'vided  into  intervals  and  the  calcu¬ 
lated  points  will  be  at  the  midpoints  of  the  intervals.  The  first 
point  will  be  1/2  interval  past  (X,  Y,  Z)ul  and  the  final  point  will 
be  within  1/2  interval  of  (X,  Y,  Z)lR.  If  (X,  Y,  Z)Lr  does  not 
lie  on  the  U-V  plane,  or  if  the  U  and  V  axes  are  not  orthogonal, 
the  location  of  the  final  point  is  not  readily  predictable. 

The  simplest  method  to  obtain  the  correct  results  is  to  specify  two 
diagonal  corners  of  the  plane  of  the  slice  on  Card  PD,  with  the  top 
having  the  short  dimension,  if  it  is  not  square,  and  the  side  having 
the  long  dimension.  Then,  on  Card  PE,  specify  the  U  axis  to  be 
parallel  to  the  edge  of  the  slice  with  the  large  dimension  (left  or 
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right  side),  and  the  V  axis  to  be  parallel  to  the  edge  of  the  slice 
with  the  small  dimension  (top  or  bottom).  Finally,  let  the  only 
entry  on  Card  PF  be  NV  equal  to  maximum  number  of  characters 
per  line  on  your  printer,  this  will  provide  the  largest  undistorted 
picture. 
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